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Abstract 

The  concept  of  Measurements  Quality  Objectives  (MQOs),  in  bioassessment 
programs  is  a  useful  tool  in  evaluating  the  consistency  of  data  and  limiting  variability  and 
potential  sources  of  measurement  error.  Typical  evaluations  of  data  repeatability 
and/or  data  quality  center  on  the  use  of  a  series  of  calculations  that  quantify  variability 
between  measures.  These  calculations  provided  some  indication  of  not  only  the  quality  of 
the  data  collected,  but  also  acted  as  a  measure  of  how  representative  the  biological  data 
were  to  each  ecoregion.  The  evaluation  of  the  Quality  Control  data  for  this  project 
provides  a  framework  for  data  users  and  water  resource  managers  to  assess  the  reliability 
and  inherent  variability  of  the  proposed  biotic  indices  for  the  state  of  Georgia. 

In  bioassessment  programs,  it  is  important  to  identify  natural  variability  of 
reference  and  impaired  sites,  as  well  as  the  variability  of  the  influences  anthropogenic 
stressors.  Calculations  of  variance  within  the  biological  parameters  measured  are 
necessary  for  identifying  the  effects  of  measurement  errors  and/or  inherent  differences 
between  sampling  sites  in  relation  to  the  overall  variance  of  a  metric  or  index  on  an 
ecoregional  and  sub-ecoregional  level. 

Considering  the  invertebrate  data  produced  by  the  Georgia  Ecoregions  Project, 
the  consistency  of  all  metric  categories  having  average  precision  measures  above  the 
prescribed  MQOs  for  both  raw  metric  values  and  standardized  metric  scores  may 
demonstrate  that  the  lotic  systems  across  the  state  of  Georgia  naturally  have  high 
variability  from  year-to-year  and  spatially  within  catchments.  This  in  turn  may  indicate 
that  the  established  precision  thresholds  of  the  MQOs  may  not  be  indicative  of  the  data 
quality  for  this  specific  project. 
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Introduction 

In  recognizing  the  need  to  improve  water  quality  conditions  of  the  surface  waters 
of  the  United  States,  mandates  have  been  set  forth  by  congress  through  the  Federal  Water 
Pollution  Control  Act,  and  more  specifically,  through  the  requirements  of  the  Clean 
Water  Act  (CWA)  (33  U.S.C.  §  et  seq.).  The  primary  objective  of  Section  101(a)  of  the 
Clean  Water  Act  is  to  "restore  and  maintain  the  chemical,  physical,  and  biological 
integrity  of  the  nation's  waters,"  (CWA,  §  101(a),  33  U.S.C,  1251(a),  1999).  Recent 
recommendations  from  the  Environmental  Protection  Agency  (EPA)  (USEPA  1987)  have 
determined  a  need  to  identify  and  establish  biological  standards  for  surface  waters.  To 
meet  the  primary  objective  of  Section  101(a),  the  Clean  Water  Act  requires  states  to 
develop  water  quality  criteria  based  on  biological  assessment  ("bioassessment").  These 
mandated  biological  criteria  ("biocriteria")  are  to  be  used  to  enforce  water  quality 
parameters  and  to  assess  possible  nonpoint  sources  (NPS)  of  pollution  (CWA,  §  319,  33 
U.S.C,  §1329,  1999). 

Biological  monitoring  has  been  mandated  by  Section  319  of  the  Clean  Water  Act 
as  the  appropriate  tool  for  assessing  the  ecological  integrity  of  streams  and  rivers 
throughout  the  United  States.  "Biological  monitoring  can  be  defined  as  the  systematic 
use  of  biological  responses  to  evaluate  changes  in  the  environment  with  the  intent  to  use 
this  information  in  a  quality  control  program.  These  changes  are  often  due  to 
anthropogenic  sources...."  (Matthews  et  at  1982)  Anthropogenic  influences  in  an 
aquatic  ecosystem  can  take  on  many  forms.  Most  obviously,  effluent  discharges  typical 
of  industry  and/or  wastewater  treatment  facilities  are  known  and  common  sources  of 
"point-source"  pollution  in  aquatic  ecosystems.  But  the  core  of  biomonitoring,  through 
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the  use  of  bioassessment  protocols,  changes  the  focus  to  NPS  pollutants  and  the 
interaction  of  widespread  chemical  or  physical  degradation  and  typical  land  use  patterns 
on  a  regional  basis.  The  effects  of  silviculture,  agriculture,  and  urban  development  on 
aquatic  resources  can  be  difficult  to  quantify,  but  through  established  biological 
monitoring  and  assessment  protocols,  the  effects  of  NPS  pollutants  can  be  identified 
(Barbour  et  al.  1996). 

In  response  to  the  bioassessment  requirements  of  the  Clean  Water  Act,  the  EPA 
has  published  two  major  guidelines  for  the  development  and  use  of  bioassessment 
protocols,  as  well  as  the  interpretation  of  the  resulting  biocriteria:  Biological  criteria: 
Technical  Guidance  for  Streams  and  Small  Rivers  (Gibson  et  al.  1996)  and  Rapid 
Bioassessment  Protocols  for  Use  in  Streams  and  Wadeable  Rivers:  Periphyton,  Benthic 
Macroinvertebrates,  and  Fish  (Barbour  et  al.  1999). 

The  Rapid  Bioassessment  Protocol  (RBP)  is  a  stepwise  methodology  for 
collecting  and  analyzing  biological,  chemical,  and  physical  habitat  data  from  stream 
ecosystems  in  order  to  provide  a  biological  framework  for  water  resource  managers  when 
assessing  water  quality  issues.  Typically,  the  result  of  a  prescribed  bioassessment 
protocol  is  the  characterization  of  "biccntena^  that  quantify  a  level  of  impairment, 
whether  minimal  or  extreme,  in  an  aquatic  system  (Fore  et  al.  1996).  Bioassessment 
protocols  and  the  resulting  biocriteria  are  an  effective  way  to  assess  water  quality  because 
of  the  integration  of  chemical  and  physical  parameters  affecting  the  biological 
community  (Karr  1990). 

The  concept  of  biocriteria  has  been  developed  to  address  the  needs  of  developing 
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biological  sampling  as  a  prime  component  of  surface  water  management  programs. 
Biocriteria  are  a  series  of  numeric  values  derived  from  the  presence  and/or  absence  of 
taxonomic  groups  in  an  ecosystem  or  region  that,  in  turn,  describe  the  biological  status  of 
the  chemical  and/or  physical  conditions  present.  Accurate  characterization  of  the 
biological  condition  involves  a  method  that  can  evaluate  patterns  of  biotic  responses  from 
the  individual  organism  to  the  ecosystem  level  (Karr  et  al.  1986). 

Traditionally,  biomonitoring  had  been  utilized  to  quantify  "before-and-after,, 
impacts  from  a  known  disturbance.  Biomonitoring,  as  currently  implemented,  can  be 
used  to  predict  impacts  to  an  aquatic  system  prior  to  major  anthropogenic  impairments 
within  a  watershed  (Rosenberg  and  Snow  1977),  as  well  as  to  serve  as  a  template, 
ensuring  compliance  to  statutory  requirements  as  set  by  the  EPA  through  the  Clean  Water 
Act.  For  compliance  measures,  biological  criteria  can  be  applied  to  evaluate  the  effects 
of  effluent  discharges  or  other  human-induced  changes  within  a  catchment  and  to 
document  that  water  quality  standards  have,  or  have  not,  been  violated  (Roper  1985). 

To  date,  the  composition  of  benthic  macroinvertebrate  communities  is  the  basis  of 
approximately  90%  of  rapid  bioassessment  programs  for  running  waters  in  the  United 
States  (Southerland  and  Stribling  1995).  The  use  of  macroinvertebrates  for 
characterizing  the  biological  integrity  of  a  waterbody  is  advantageous  for  a  number  of 
reasons.  The  large  numbers  of  invertebrate  species  present  in  aquatic  systems  provide  an 
assortment  of  biological  responses  to  induced  stresses  (Hellawell  1986).  The  sessile 
nature  and  limited  dispersal  ability  of  aquatic  macroinvertebrates  make  them  ideal 
indicators  of  not  only  ambient  water  quality  conditions  (Hawkes  1979),  but  also  as  a 
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gauge  for  spatial  and  temporal  environmental  disturbances,  both  natural  and 
anthropogenic. 
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The  physical  and  chemical  changes  that  can  occur  in  a  catchment,  as  a  result  of 
anthropogenic  influences,  will  have  direct  and  indirect  effects  not  only  on  the  aquatic 
communities  present,  but  also  on  the  habitat  structure  and  the  functional  food  web  of  the 
lotic  ecosystem  (Townsend  and  Riley  1999).  Although  some  benthic  macroinvertebrate 
taxa  are  widespread  throughout  the  stream  ecosystem  and  utilize  many  different  habitats, 
there  are  also  many  groups  that  are  more  restricted  to  one  specific  habitat  (Pardo  and 
Armitage  1997).  Having  distinctive  habitat  and  biological  requirements, 
macroinvertebrate  assemblages  can  be  used  in  a  predictive  manner  in  biomonitoring 
programs.  The  absence  and/or  presence  of  certain  taxonomic  groups  and  species  within  a 
stream  can  be  indicative  of  pollution  levels  in  both  impaired  and  unimpaired  catchments 
(Ravera  2001;  Cairns  and  Pratt  1993). 

Through  the  use  of  bioassessment  protocols,  an  array  of  biological  metrics  are 
developed  to  characterize  typical  aquatic  communities,  representing  both  minimally 
impaired  ('"reference")  streams  and  impaired  streams.  The  "multi-metric"  approach  to 
characterizing  water  quality  is  effective  because  of  its  analysis  of  a  number  of  biological 
responses  in  macroinvertebrate  communities.  However,  there  are  a  number  of  factors  to 
consider  when  identifying  those  metrics  that  are  most  characteristic  of  the  biological 
condition  in  stream,  catchment,  or  region. 

The  most  difficult  part  of  developing  a  multi-metric  approach  for  assessing  water 
quality  parameters  for  a  stream,  catchment,  and/or  ecoregion,  is  to  determine  which 
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biological  metrics  are  diagnostic  of  the  responses  of  benthic  macroinvertebrate 
assemblages  to  anthropogenic  changes  within  that  catchment.  The  initial  approach  to 
resolving  this  problem  requires  a  delineation  of  biological,  chemical,  and  morphological 
variability  across  a  geographic  area.  Omernik  (1987)  has  developed  an  aquatic  and 
terrestrial  map  of  ecoregions  for  the  United  States.  This  map  provides  the  framework  for 
grouping  ecosystems,  based  upon  patterns  of  topography,  geology  and  soil,  and  land  use. 
This  ecoregional  grouping  is  intended  to  minimize  variability  within  similar  regions,  as 
well  as  maximize  variability  between  dissimilar  regions. 

The  theory  behind  the  use  of  the  ecoregion  concept  is  that  adjoining  land  forms 
with  similar  geologic  features,  soil  types,  vegetation,  and  climatic  influences  will  most 
likely  possess  similar  biological  communities  (Omernik  1995;  Hughes  1995;  Omernik 
and  Gallant  1990).  This  concept  is  useful  in  conjunction  with  bioassessment  programs 
because  it  can  be  used  to  characterize  and  predict  natural  variations  among  systems 
within  similar  geographic  regions,  as  well  as  to  detect  responses  to  disturbances  based  on 
some  reference  condition  (Hughes  and  Larsen  1988). 

With  the  variable  geology  and  vegetation  patterns  across  the  state  of  Georgia,  it 
should  be  expected  that  a  variety  of  macroinvertebrate  assemblages  will  reflect  the 
ambient  water  quality  and  habitat  structure  of  those  systems.  Similarly,  any  degradation 
of  habitat  and  deviation  from  typical  water  quality  in  a  region  should  be  reflected  by 
changes  in  the  composition  of  the  macroinvertebrate  community.  Characterizing  a 
representative  macroinvertebrate  community  in  minimally  impaired  catchments  serves  as 
a  reference  point  for  other  stream  ecosystems  that  have  been  subjected  to  some  sort 
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anthropogenic  stress. 

A  reference  condition,  as  prescribed  by  bioassessment  protocols,  is  defined  as 
"the  condition  that  is  representative  of  a  group  of  minimally  disturbed  sites  organized  by 
physical,  chemical  and  biological  characteristics"  (Reynoldson  et  al.  1997).  The 
biological  condition  of  a  stream,  or  group  of  streams,  that  are  classified  as  "reference" 
then  serve  as  the  point  of  comparison  for  all  other  streams  within  a  catchment  and/or 
ecoregion.  The  chemical,  physical,  and  biological  attributes  of  a  reference  stream  can 
then  be  used  to  identify  levels  of  impairment  in  streams  that  are  known  to  be  altered.  The 
differences  between  the  biological  condition  of  a  reference  and  impaired  site  can  be 
quantified  through  a  series  of  biological  metrics.  These  metrics  are  then  used  to  develop 
a  ranking  system  to  identify  streams  that  have  acceptable  or  degraded  water  quality  per 
EPA  standards. 

To  accurately  assess  the  effects  of  anthropogenic  influences,  natural  variability 
within  these  geographical  boundaries  must  be  characterized.  In  Georgia,  there  is  a  very 
distinctive  geological,  vegetative,  and  geomorphological  transition  from  the  northwest 
region  to  the  southeast  region  {see  Table  1).  This  change  in  ecoregional  character 
dictates  a  variety  of  stream  morphologies  with  variable  habitat  structures  and  water 
chemistries.  The  final  determination  of  a  series  of  metrics  must  somehow  account  for 
natural  biological  variability  within  and  across  ecoregional  boundaries  (MDEQ  2003). 

In  1999,  Columbus  State  University  (CSU)  was  selected  to  create  a  rapid 
bioassessment  program  for  the  state  of  Georgia.  Funding  for  this  project  was  provided 
through  a  grant  from  the  Georgia  Department  of  Natural  Resources,  Environmental 
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Table  1.  Descriptions 

of  the  primary 

Georgia  ecoregions 

Ecoregion  Code  and 

Principal  Land  Use/ 

Ecoregion  Name 

Geology 

Climate 

Dominant  Vegetation 

45 

Piedmont 

Metamorphic 

Mesic-Xeric 

Silviculture  and  Urban 
Mixed  Forest 

65 

Southeastern  Plains 

Sedimentary 

Mesic-Xeric 

Agriculture  and  Silviculture 

(Cretaceous- 

Pine  Forest 

Miocene) 

66 

Blue  Ridge 
Mountains 

Metamorphic 

Mesic-Submesic 

Hardwood  Forest 

67 

Ridge  and  Valley 

Sedimentary 

Mesic-Submesic 

Agriculture 

(Paleozoic) 

Hardwood  Forest 

68 

Cumberland  Plateau 

Sedimentary 

Mesic-Submesic 

Agriculture 

(Paleozoic) 

Hardwood  Forest 

75 

Southern  Coastal 
Plain 

Sedimentary 

Mari-time 

Agriculture  and  Silviculture 

(Pliocene- 

Pine  Forest 

Pleistocene) 

Ecoregion  delineation  per  Omernik  (1987)  and  ecoregion  descriptions  per  Wharton 
(1989). 


Protection  Division  (GAEPD),  via  sponsorship  of  the  United  States  Environmental 
Protection  Agency  (more  specifically,  United  States  Environmental  Protection  Agency 
Clean  Water  Act,  Section  319(h)  FY  98  -  Element  1)  funding.  The  resulting  "Georgia 
Ecoregions  Project"  consisted  of  four  phases  of  biological,  chemical,  and  physical 
data    collection  to  characterize    water  quality  conditions  across  the  state.    The  final 
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analysis  of  this  data  will  be  used  to  establish  biocriteria  relevant  to  each  geographic 
region  of  Georgia. 

Ultimately,  the  resulting  biocriteria  and  numerical  ranking  system  derived  from 
this  project  can  be  used  to  evaluate  the  possible  sources  and  effects  of  NPS  pollution  and 
the  effectiveness  of  Best  Management  Practices  (BMPs)  to  control  NPS  inputs,  as  well  as 
to  assess  the  level  of  impacts  from  Total  Maximum  Daily  Loads  (TMDLs)  in  aquatic 
systems.  These  established  biocriteria  can  also  be  used  to  identify  regions  or  specific 
catchments  in  need  of  restoration,  as  well  as  to  characterize  the  sources  of  impairment 
and  to  monitor  trends  over  time  (Barbour  et  al.  1999). 

Through  the  use  of  Geographical  Information  Systems  (G1S)  software,  {i.e. 
Arc  View),  land  use  data,  and  best  professional  judgment,  efforts  were  made  to  identify 
and  locate  as  many  potential  "reference1'  and  "impaired"  sampling  sites  as  possible 
(Olson  2002,  Gore  et  al.  2004).  This  process  was  necessary  in  order  to  adequately 
illustrate  the  inherent  biological,  chemical,  and  physical  variability  of  streams  throughout 
the  state  within  their  catchments  and  ecoregional  boundaries.  The  goal  was  to  collect 
physical  and  chemical  data  from  a  minimum  of  ten  streams  specific  to  each  sub- 
ecoregion  identified  for  Georgia,  five  sites  being  classified  as  reference,  (or  minimally 
impaired),  and  five  sites  being  classified  as  impaired,  based  on  land  use  parameters 
within  the  catchment.  A  statistical  summary  of  land  use  within  the  ecoregions  of  Georgia 
is  provided  in  Appendix  A. 

In  conjunction  with  these  proposed  reference  and  impaired  sites  for  sampling, 
additional  samples  were  collected  as  dictated  by  the  Quality  Assurance  Project  Plan 
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(QAPP)  (CSU  2000).  Throughout  each  phase  of  the  Georgia  Ecoregion  Study,  there 
were  a  number  of  duplicate  samples  taken  to  satisfy  the  Quality  Assurance/Quality 
Control  (QA/QC)  requirements  of  the  QAPP.  These  duplicate  samples  were  taken  in 
order  to  assess  the  repeatability  and  precision  of  the  collected  data,  as  well  as  to  assess 
the  training  and  level  of  effort  between  and  among  field  teams.  In  this  paper,  QC  data  are 
assessed  in  terms  of  Measurement  Quality  Objectives  (MQOs)  as  outlined  by  the  QAPP, 
but  address,  more  specifically,  the  amount  and  degree  of  variability  present  in  the 
wadeable  stream  ecosystems  across  the  state  of  Georgia,  and  how  these  samples  affect 
initial  characterization  of  the  reference  condition. 

There  were  two  designations  for  QC  samples  collected  to  satisfy  the  QAPP 
document:  (1)  a  spatial,  "duplicate  reach"  QC  sample,  and  (2)  a  temporal,  "phase"  QC 
sample.  According  to  QAPP  procedures,  the  QC  type  of  duplicate  sampling  is  performed 
in  order  to  assess  the  precision  and  accuracy  of  the  field  teams  and  the  representativeness 
of  the  data  as  some  measure  of  "data  quality"  in  bioassessment  programs.  It  is  important 
to  analyze  the  consistency  of  field  teams  to  ensure  that  personnel  are  properly  trained  so 
that  the  collection  of  biological  data  are  free  from  bias  and  error,  but  more  importantly, 
for  this  paper,  the  objective  of  analyzing  the  additional  biological  data  was  to  determine 
if  the  sites  chosen  to  characterize  the  biological  condition  were  true  representations  of  the 
biological  community  in  that  stream. 

During  development  of  the  reference  condition  for  Georgia,  the  additional  data 
collected  through  the  QC  samples  were  not  used  in  the  creation  of  overall  metric  scores, 
or  in    characterization  of  the    final    biological  index.    Thus,  it  became  important  to 
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determine  the  point  that  the  data  thoroughly  and  accurately  reflected  the  composition  of 
the  macroinvertebrate  community.  There  have  been  numerous  studies  of  the  effect  of 
sample  size  on  the  variability  of  biotic  indices  in  bioassessment  programs  (Li  et  al.  2001; 
Metzeling  and  Miller  2001).  Increase  in  sample  size  will  result  in  the  increase  of  number 
of  individuals  collected,  but,  more  importantly,  also  corresponds  to  an  increase  in  the 
number  of  taxa  in  the  system  being  sampled.  It  has  also  been  demonstrated  that 
increasing  the  size  of  the  sampling  area,  (whether  it  be  sampling  more  than  one  riffle  or  a 
combination  of  habitats  to  constitute  one  sample),  has  an  effect  on  the  range  of  variance 
(Hannaford  and  Resh  1995,  Norris  et  al.  1993).  When  considering  bioassessment 
protocols,  it  raises  the  question  of  determining  what  important  taxa  may  have  been  be 
missed  and  how  these  excluded  taxa  may  influence  the  range  of  variability  of  the  metrics 
used  to  determine  the  reference  and  impaired  condition. 

There  have  been  a  number  of  papers  analyzing  variability  in  data  using  RBP 
protocols  (see  reviews  by  Hannaford  and  Resh  1995),  but  the  majority  of  these  have 
centered  on  specific  habitat  types  such  as  riffles  and  runs  (see,  for  example,  Feminella 
2000),  based  upon  the  assumption  that  swifter  water  habitats  yield  the  highest  species 
richness  and  abundance  of  invertebrates  (Hynes  1970;  Allan  1995).  Also  common  in 
these  previous  studies  has  been  the  use  of  "in-field"  subsampling  of  macroinvertebrates 
as  the  basis  for  characterizing  variability  in  the  data  sets  (Metzeling  and  Miller  2001). 
Logically  there  is  some  question  about  bias  resulting  from  "in-field"  subsampling  of 
macroinvertebrates,  as  there  may  be  a  tendency  to  choose  the  larger,  more  obvious 
organisms  for  analysis,  resulting  in  skewed  final  metrics  and  biotic  indices  calculated  for 
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a  stream. 

Additionally,  given  that  many  macroinvertebrates  have  very  specific  habitat 
requirements,  it  is  to  be  expected  that  metric  results  would  vary  as  a  function  of  the  range 
of  particular  habitats  being  sampled.  There  are  numerous  species  that  thrive  in  habitats 
such  as  tree  roots  along  stream  banks  and  woody  debris  (i.e.  snags),  an  especially 
productive  habitat  type  in  low-gradient  stream  systems  typical  of  southern  Georgia 
(Benke  et  al.  1985).  The  sampling  of  multiple  habitats  in  bioassessment  protocols 
provides  a  better  biological  "picture'"  of  the  faunal  communities  that  are  subject  to 
changes  in  habitat  structure  and  water  quality. 

The  research  described  here  addressed  a  number  of  questions.  With  regards  to  the 
Georgia  Ecoregions  QC  data,  does  the  inclusion  of  additional  taxonomic  data  change  the 
range  of  variability,  and  what  criteria  define  the  reference  or  impaired  condition? 
Second,  will  the  restriction  or  expansion  of  those  ranges  of  variability  create  difficulties 
in  interpretation  of  anthropogenic  stressors  on  the  biotic  community?  Likewise,  does  the 
range  of  variability  within  the  identified  metrics  and  biotic  indices  hinder  the  decision 
making  process  for  water  resource  managers?  The  answers  to  these  questions  might 
indicate  that  increasing  the  sample  size,  (e.g.  increasing  the  number  of  reference  and 
impaired  sites  samples,  and/or  increasing  the  reach  length),  in  RPB  programs  may  better 
characterize  natural  variability  within  and  between  ecoregions,  and  also  reduce  the 
variability  of  the  final  metrics  used  to  characterize  the  reference  condition  and  water 
quality,  as  well  as  more  narrowly  defining  numerical  criteria  of  stream  health  (Gore  et  al. 
2005). 


.       ■ 
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Materials  and  Methods 

The  Georgia  Ecoregions  Project  consisted  of  four  phases  of  ecoregional 
delineation,  sample  collection,  and  data  analysis.  Phase  1  consisted  of  identifying  and 
delineating  ecoregional  boundaries,  as  established  by  Omernik  (1987).  The  sub- 
ecoregions  used  in  this  study  were  from  the  Level  III  and  IV  Sub-Ecoregions  of  Georgia 
(Griffith  et  al.  2001)  as  illustrated  in  Figure  1.  Phase  2  consisted  of  evaluating  candidate 
reference  streams  based  upon  abiotic  factors  (i.e.  surrounding  land  use  patterns,  physical 
habitat  quality,  and  water  chemistry).  In  that  phase,  criteria  were  established  to 
characterize  a  "reference'''  stream  in  terms  of  water  quality  and  biological  (primarily, 
macroinvertebrate)  assemblages  (Olson  2002).  The  reference  condition  was  essential  to 
provide  a  "'benchmark"  by  which  impairment  status  would  be  characterized. 

In  Phases  3  and  4  of  this  project,  the  process  of  identifying  and  sampling 
reference  and  impaired  sites  was  continued  in  order  to  collect  as  much  biological  data  as 
possible  and  accurately  reflect  water  quality  and  macroinvertebrate  assemblages 
characteristic  of  the  defined  ecoregions.  To  minimize  the  effect  of  temporal  variability 
between  sampling  years/phases,  a  sampling  season  (or  "index  period")  between  August 
and  February  was  determined  to  be  the  most  indicative  of  the  aquatic  communities  for  the 
Georgia  Ecoregions  Project. 

At  each  sampling  site,  a  series  of  physical  and  chemical  sampling  protocols  were 
performed  to  collect  data  relevant  to  the  habitat  quality,  water  quality,  and 
macroinvertebrate  communities  that  are  representative  of  the  biological  condition  of  the 
site,  as  well  as  being  indicative  of  ecoregional  character,  (i.e.  geology,    vegetation, 
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Figure  1.  Level  III  and  IV  Sub-Ecoregions  of  Georgia. 
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climate  etc.) ,  and  land  use  patterns,  {i.e.  silviculture,  urbanization,  etc.).  Chemical  and 
biological  sampling  of  the  selected  sites  was  preformed  using  the  following  procedures. 
These  sampling  protocols  are  described  further,  in  more  detail  in  Columbus  State 
University's  (CSU)  QAPP  as  Standard  Operating  Procedures  (SOPs)  and  generally 
follow  the  recommendations  of  the  RBP  (Barbour  et  al.  1999): 

1.  Benthic  macroinvertebrate  communities  were  collected  using  the  Georgia  DNR's 
20-Jab  Method  (CSU  2000).  Table  2  summarizes  the  level  of  effort  prescribed  for 
various  habitat  types  that  are  characteristic  of  high-  and  low-gradient  stream 
systems.  Macroinvertebrates  collected  from  these  habitat  types  were  composited 
into  a  single  sample  and  returned  to  the  lab  for  further  processing. 
Macroinvertebrates  were  identified  to  the  lowest  practicable  taxonomic  level  and 
enumerated. 

2.  Water  chemistry  was  measured  both  in  situ,  using  a  Hydrolab  H-20  probe,  and  by 
grab  samples  that  were  analyzed  at  a  later  time  in  the  CSU  laboratory.  The  water 
chemistry  parameters  that  were  analyzed  for  this  project  are  listed  in  Appendix  B. 

3.  The  physical  properties  of  the  streams  were  also  recorded.  Those  properties 
included  a  streambed  cross  section,  velocity,  substrate  size  and  composition  using  a 
modified  Wolman  Pebble  Count  (Bevenger  and  King  1995).  Additional 
observations  included  the  extent  of  canopy  cover,  presence  of  oils  and/or  odors, 
adjacent  land  use  along  the  stream  channel,  bank  erosion,  and  types  of  deposits,  {see 
Appendix  C). 
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4.  Visually-based  habitat  assessments  were  also  completed  for  each  site  using  the 
EPA's  Rapid  Bioassessment  Protocol  habitat  assessment  methods  and  forms,  (see 
Appendix  D).  Similar  to  the  partitioning  of  habitat  types  between  stream  gradients 
for  the  sampling  of  macroinvertebrate  communities,  the  habitat  assessment  forms 
used  to  characterize  physical  and  geomorphological  attributes  of  a  stream  ecosystem 
were  also  categorized  by  gradient  classification. 


Table  2*.    Prioritized  list  of  habitat  types  for  sampling  and  « 
modified  20-jab  method. 

sample  reallocation  for  the 

HIGH  GRADIENT  STREAMS 

Priority 

Habitat  Type 

Number  of  Samples 

1 

Fast  Riffle 

3 

2 

Slow  Riffle 

3 

3 

Snags 

5 

4 

Undercut  Banks/Rootwads 

3 

5 

Leaf  Packs 

3 

6 

Sand 

3 

7 

Macrophytes  (if  any) 

3 

LOW  GRADIENT  STREAMS 

Priority 

Habitat  Type 

Number  of  Samples 

1 

Woody  debris/Snags 

8 

2 

Undercut  Banks/Rootwads 

6 

3 

Leaf  Packs 

3 

4 

Sand 

3 

5 

Macrophytes  (if  any) 

3 

*From    Columbus    State    University's    Quality    Assurance    Project    Plan    document 
(CSU2000). 


Each  macroinvertebrate  sample  collected  in  the  field  was  preserved  in  70% 
ethanol  until  further  processing.  The  macroinvertebrate  samples  were  evenly  spread  out 
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on  a  Caton  gridded  screen  for  subsampling  (Caton  1991).  Using  random  number 
generation,  squares  within  the  grid  system  were  designated  for  removal  and  subsequent 
examination.  Each  grid  square  was  placed  in  a  "white-pan"  and  examined  for  the 
presence  of  macroinvertebrates.  All  organisms  from  the  square  were  removed  and  again 
preserved  in  70%  ethanol.  For  each  macroinvertebrate  sample,  the  goal  was  to  collect  a 
total  of  200  organisms.  Squares  for  the  grid  were  continuously  selected  until  the  required 
number  of  organisms  was  collected.  Macroinvertebrates  were  then  identified  to  the 
lowest  possible  taxonomic  level  and  entered  into  a  database  developed  by  Tetra  Tech 
(1999)  known  as  the  Ecological  Data  Application  System  (EDAS)  for  metric  analysis. 

The  QAPP  describes  the  procedures  that  were  used  in  data  collection  and  their 
rationale,  as  well  as  a  series  of  activities  and  reporting  procedures  that  were  used  to 
document  data  quality.  As  prescribed  by  the  QAPP  document  for  the  Georgia  Ecoregions 
Project,  a  number  of  sites  were  designated  for  additional  sampling.  To  address  QC/QA 
protocols  related  to  data  quality,  ten  percent  of  all  the  designated  sampling  sites  were 
required  to  have  duplicate  sampling  performed.  These  duplicate  samples  fell  into  two 
designations:  "spatial"  (200  meter  QC)  and  "temporal"  (Phase  QC). 

Quality  control  samples  that  were  designated  as  "spatial"  essentially  "doubled" 
the  length  of  the  reach  designated  for  sampling.  Once  the  primary  sample  reach  of  one 
hundred  meters  was  established,  and  all  RBP  sampling  requirements  satisfied,  the 
immediate,  next  one  hundred  meter  reach  was  sampled.  "Temporal"  QC's  were  sites  that 
were  sampled  in  succeeding  phases  of  the  ecoregions  project,  where  the  originally 
established  sample  site    reach  was  resampled  in  a  subsequent    "index  period".  This 
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sampling  approach  addressed  two  possible  variations  within  a  stream  ecosystem:  (1)  the 
variability  of  the  distribution  of  habitats  longitudinally  within  a  catchment,  and  (2) 
changes  of  the  macroinvertebrate  communities  over  time. 

All  QC  sites  were  randomly  chosen  via  a  random  number  generation  function  in 
Microsoft's  Excel  program.  As  a  result  of  this  random  number  generation,  there  was 
some  unevenness  in  the  number  of  duplicate  reference-  and  impaired-site  samples,  as 
well  as  the  number  of  spatial-  and  temporal-QC  samples  collected.  Additionally,  the  total 
number  of  QC  sites  collected  for  this  project  was  not  evenly  distributed  throughout  each 
ecoregion  and  subecoregion.  A  list  of  all  spatial  and  temporal  QCs  samples  collected  per 
ecoregion  is  provided  in  Appendix  E,  combined  with  GIS  maps  to  illustrate  their 
geographic  locations  throughout  the  state  of  Georgia. 

As  sites  were  sampled  and  taxonomic  identifications  were  completed,  all  of  the 
physical,  chemical,  and  biological  data  gathered  were  entered  into  the  EDAS  database  for 
further  analysis.  The  calculated  macroinvertebrate  metrics  encompass  a  number  of 
benthic  macroinvertebrate  community  structures  and  functions  that  characterize  the 
ecological  status  of  the  aquatic  system  being  analyzed.  An  assortment  of  approximately 
65  metrics,  from  five  major  structural  and  functional  groups  (i.e.,  taxonomic  richness, 
community  composition,  tolerant/intolerant  organisms,  functional  feeding  groups,  and 
life  habit)  were  calculated  for  each  stream  sampled.  A  list  of  all  metrics  considered  in 
developing  biocriteria  for  the  Georgia  Ecoregions  Project  has  been  compiled  in  Table  3. 
A  brief  description  of  each  metric  groups  and  its  significance  to  characterizing  ambient 
water  quality  conditions  follows: 


• 


• 


18 

•  Taxonomic  Richness  -  metrics  included  in  this  group  evaluated  the  number  of 
individual  taxa  within  larger  taxonomic  groups  (i.e.  number  of  families  or  genera 
within  an  order  of  aquatic  invertebrates  such  as  Ephemeroptera,  Diptera,  etc. ). 
Typically,  high  values  of  taxonomic  richness  are  indicative  of  better  water  quality 
and  a  healthier  lotic  ecosystem. 

•  Community  Composition  -  these  metric  values  are  expressed  as  percentages, 
representing  a  proportion  of  individuals  in  a  sample  belonging  to  some  specific 
taxonomic  group.  Higher  percentages  of  those  organisms  that  have  been  known 
to  tolerate  degraded  conditions  (i.e.  Diptera)  are  assumed  to  be  characteristic  of 
impaired  water  quality. 

•  Tolerant/Intolerant  Taxa  -  this  group  of  metrics  are  represented  by  the 
tolerance  levels  of  biota  to  stress.  In  systems  with  high  anthropogenic  stress,  taxa 
classified  as  intolerant  to  pollution  are  assumed  to  be  the  first  organisms  to  be 
eliminated  from  the  ecosystem,  becoming  less  abundant.  Concurrently,  those  taxa 
with  higher  tolerance  to  pollution  impacts  are  assumed  to  dominate  the  system. 

•  Functional  Feeding  Group  -  these  metrics  reflected  the  dominant  feeding  mode 
of  the  biological  community  in  the  sample.  The  ecological  responses  of 
organisms  with  specialized  or  generalized  feeding  habits  are  assumed  to  be 
indicative  of  pollution  or  anthropogenic  disturbances.  For  example,  the 
abundance  of  "shredders"  and  "filterers"  can  be  affected  when  organic  materials 
become  scarce  or  unsuitable. 
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Life  Habit  -  the  metrics  include  measures  of  taxa  richness  and  composition  that 
described  the  locomotive  and  positioning  mechanisms  of  benthic 
macroinvertebrates  (i.e.  burrowing,  swimming,  etc.).  This  group  of  metrics  are 
probably  the  most  difficult  to  characterize  responses  to  anthropogenic  stressors,  as 
there  has  been  no  definitive  indication  of  how  these  communities  identified  by  life 
habit  respond  to  increased  or  decreased  perturbations  in  the  ecosystem.  (Kerans 
and  Karr  1994). 


Table  3.  Predicted  Responses  of  Benthic  Macroinvertebrate  Metrics  to  Stress. 


METRIC  CATEGORY 

METRIC 

STRESS 
RESPONSE 

Taxonomic  Richness 

Total  Taxa 

Decrease 

Ephemeroptera, 

Plecoptera,  & 

Trichoptera  (EPT)  Taxa 

Decrease 

Ephemeroptera  Taxa 

Decrease 

Plecoptera  Taxa 

Decrease 

Trichoptera  Taxa 

Decrease 

Coleoptera  Taxa 

Decrease 

Diptera  Taxa 

Decrease 

Chironomidae  Taxa 

Decrease 

Tanytarsini  Taxa 

Decrease 

Evenness 

Decrease 

Margalef  s  Index 

Decrease 

Shannon- Wiener  base  e 

Decrease 

Simpson's  'Diversity 

Increase 

METRIC  CATEGORY 

METRIC 

STRESS 
RESPONSE 

Community  Composition 

%  Ephemeroptera 

Decrease 

%  Amphipoda 

Decrease 

%  Chironomidae 

Increase 

%  Coleoptera 

Decrease 

%  Diptera 

Increase 

%  Gastropoda 

Decrease 

. 
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Table  3.     Predicted  Responses  of  Benthic  Macroinvertebrate  Metrics  to  Stress, 
(cont) 

METRIC  CATEGORY 

METRIC 

STRESS 
RESPONSE 

Community  Composition 
(cont.) 

%  Isopoda 

Increase 

%  Nonlnsect 

Increase 

%  Odonata 

Increase 

%  Plecoptera 

Decrease 

%  Tanytarsini 

Decrease 

%  Oligochaeta 

Increase 

%  Trichoptera 

Decrease 

%  Chironominae  /  Total 
Chironomidae  (TC) 

Variable 

%  Orthocladiinae  /  TC 

Decrease 

%  Tanypodinae  /  TC 

Increase 

%  Hydropsychidae  / 
Total  Trichoptera 

Increase 

%  Hydropsychidae  / 
Total  EPT 

Increase 

%  Tanytarsini  /  TC 

Decrease 

%  Cricotopus  sp.  & 
Chironomus  sp.  1  TC 

Increase 

METRIC  CATEGORY 

METRIC 

STRESS 
RESPONSE 

Tolerance/Intolerance 

Tolerant  Taxa 

Increase 

%  Tolerant  Individuals 

Increase 

Intolerant  Taxa 

Decrease 

%  Intolerant  Individuals 

Decrease 

%  Dominant  Individuals 

Increase 

Dominant  Individuals 

Increase 

Beck's  Index 

Decrease 

Hilsenhoff  s  Biotic  Index 
(HBI) 

Increase 

North  Carolina  Biotic 
Index  (NCBI) 

Increase 

METRIC  CATEGORY 

METRIC 

STRESS 
RESPONSE 

Functional  Feeding  Group 

%  Scraper 

Decrease 

Scraper  Taxa 

Decrease 

%  Collector 

Decrease 

Collector  Taxa 

Decrease 

%  Predator 

Decrease 

Predator  Taxa 

Decrease 
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Table  3.     Predicted  Responses  of  Benthic  Macroinvertebrate  Metrics  to  Stress, 
(cont.) 


METRIC  CATEGORY 

METRIC 

STRESS 
RESPONSE 

Functional  Feeding  Group 
{cont.) 

%  Shredder 

Decrease 

Shredder  Taxa 

Decrease 

%  Filterer 

Increase 

Filterer  Taxa 

Decrease 

METRIC  CATEGORY 

METRIC 

STRESS 
RESPONSE 

Life  Habit 

Clinger  Taxa 

Decrease 

%  Clinger 

Decrease 

Burrower  Taxa 

Decrease 

Climber  Taxa 

Decrease 

Sprawler  Taxa 

Decrease 

Swimmer  Taxa 

Decrease 

The  use  of  various  metrics  from  these  groups  resulted  in  a  "multi-metric" 
approach  to  assess  the  health  of  a  stream  ecosystem.  This  variety  of  biological  data 
includes  many  ecologically  significant  factors  in  aquatic  systems  that  are  then  compiled 
into  a  single  biotic  index  relevant  to  each  ecoregion.  Final  biotic  indices  were  comprised 
of  five  to  seven  metrics,  with  the  metrics  being  chosen  with  at  least  one  representative 
from  each  of  the  five  metric  categories  mentioned  above. 

Metric  values  calculated  by  EDAS  were  separated  by  ecoregion  and 
subecoregion,  as  well  as  by  impairment  status,  within  the  ecoregional  designation.  These 
raw  metric  scores  were  initially  used  to  distinguish  which  metrics  were  to  be  considered 
as  candidates  for  the  final  biotic  index.  A  series  of  analyses  was  performed  to  assess  the 
ability  of  each  metric  to  discriminate  between  the  characteristics  of  a  "reference"  stream 
and  an  "impaired"  stream.  For  those  sites  designated  as  reference  streams,  inter-quartile 
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ranges  of  the  distribution  of  metric  values  for  each  metric  category  were  calculated. 
These  percentile  values  then  served  as  templates  for  the  calculations  of  "discrimination 
efficiency",  metric  score  standardization,  and  the  final  biotic  index  for  each  ecoregion 
and  subecoregion. 

Raw  metrics  scores  for  both  reference  and  impaired  streams  sites  were  used  to 
calculate  the  discrimination  efficiency  (DE)  of  each  metric.  DE  was  calculated  by  the 
following  formula  (Stribling  et  al.  2000): 

DE  =  (a/b)  x  100 
For  those  metrics  values  that  typically  decrease  with  increased  stress,  a  =  the  number  of 
impaired  streams  that  scored  below  the  25th  percentile  of  the  distribution  of  reference 
metric  values;  whereas  for  those  metrics  values  that  increase  in  response  to  increasing 
stress,  a  =  the  number  of  impaired  streams  that  scored  above  the  75th  percentile  of  the 
distribution  of  the  reference  metric  values,  and  b  =  the  total  number  of  impaired  sites 
sampled  specific  to  the  ecoregion  and/or  subecoregion  analyzed.  No  metric  with  a  DE 
value  below  0.5  was  considered  as  a  candidate  for  inclusion  in  the  final  biotic  index. 

Those  metrics  that  exhibited  high  DEs  were  then  evaluated  with  box-and-whisker 
plots  (Tukey  1970)  to  compare  the  distribution  of  metric  values  for  reference  and 
impaired  stream  classes.  Box-and-whisker  plots  graphically  demonstrate  the  range  of 
variation  not  only  within  the  stream  class  {i.e.  reference  or  impaired),  but  also  between 
the  stream  classes  (i.e.  reference  vs.  impaired)  (Gibson  et  al.  1996).  The  best  performing 
metrics,  and  those  that  were  to  be  considered  as  candidates  for  the  final  index, 
demonstrated  minimal  variation  within  stream  classes  and  maximized  variation  between 
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stream  classes  (Barbour  et  al.  1999).  The  metrics  deemed  to  be  candidates  were  then 
analyzed  by  Pearson-Product-Moment-Correlation  to  determine  which  metrics  may  be 
redundant  with  one  another  (Steel  and  Torrie  1980).  Redundant  metrics  were  not  used 
when  determining  the  final  biotic  index  for  each  ecoregion. 

To  incorporate  raw  metric  values  into  an  index,  the  data  were  first  standardized, 
thus  creating  a  metric  "score".  Standardization  of  the  raw  metric  values  normalized  each 
metric  to  a  similar  scale  so  that  they  could  be  compared  to  one  another.  The  method  of 
standardization  of  the  raw  metric  values  was  dependent  upon  the  stress  response  of  the 
metric.  For  metric  values  that  were  known  to  decrease  as  a  response  increasing  stress, 
the  standardized  score  was  calculated  as: 

Standardized  Metric  Score  =  (c/d)  x  100 
Where  c  =  the  raw  metric  value  for  a  site,  and  d  =  the  95th  percentile  of  the  distribution  of 
the  reference  metric  values.  For  metric  values  that  were  known  to  increase  as  a  response 
to  increasing  stress,  the  standardized  score  was  calculated  as: 

Standardized  Metric  Score  =  |(e-c)/(e-f)]  x  100 
Where  c  =  the  raw  metric  value  for  a  site;  e  =  the  maximum  observed  value  among  all 
streams  {i.e.,  reference  and  impaired);  and  f  =  the  5th  percentile  of  the  distribution  of  the 
reference  metric  values. 

Standardized  metric  scores  of  the  initial  candidate  metrics  were  then  combined 
from  each  metric  category  to  create  a  candidate  biotic  index.  A  series  of  indices  was  first 
calculated  to  see  which  combination  of  metrics  determined  the  best  index  to  describe  the 
biotic    factors  that    represented    the    reference  and    impaired    condition  for    specific 
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ecoregions.  Index  Scores  were  assembled  from  the  standardized  scores  in  an  additivc 
manner  as: 

Index  Score  ■  (g+ri+i+j+k...)/n 
Where  tf,  h,  i,  j,  k...  the  standardized  metric  score  ol  the  candidate  metrics  and  n  the 
total  number  of  standardized  metrics  included  in  the  index  All  candidate  indit.es 
contained  live  to  seven  metrics  that  are  scored  on  a  0  to  100  point  scale  I  he  initial 
candidate  indices  were  then  analyzed  lor  1)1  s  and  compared  b>  box-and-whiskcr  plots 
The  final  candidate  index  was  distinguished  by  the  highest  1)1  and  the  Ix-si  box-and- 
whisker  separation  A  final  biotic  index  was  de\elo|K*d  lor  each  ecoregion  and 
subecoregion  of  Georgia  (Gore  el  <//.  2004,  Middleton  2005)    Although  the  final  biotic 

index  was  guided  by  statistical  considerations  {i.e.,  1)1  s  intcr-cjuart ilc  ranges,  I  U  I  the 
choices  of  the  final  metrics  were  also  based  ii|x>n  their  ability  to  accuratelv  represenl 
ecoregional  characteristics  and  the  abihtv  to  typify  a  response  to  anthropogenic  stress 
(Barbour  el  al.  1999) 

I  here  are  main  important  factors  to  consider  when  developing  biocnteria  tor 
stream  ecosystems  In  the  process  ol  collecting  biological  data,  field  methods  cannot 
predict  if  the  information  being  collected  is  an  accurate  portrayal  of  the  ecosystem  under 
investigation  (I  I  I'M  1995)  I  he  properties  ol'  a  given  field  sample  ma\  be  known,  but 
tvpicallv  biological  data  are  collected  with  the  intent  of  answering  questions  relalu 

much  largei  spatial  and  temporal  scales  (Barboui  ./  <v/.  1999)     I  he  consistency  of  Odd 

methods   and    level    of  effort    in   collecting   biological    data   aic   the    kev    to   obtaining 

information  thai  is  representative  o\  field  conditions  at  that  point  in  tune,  but  trulv 
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accurate  assessments  of  the  biological  data  are  hindered  because  the  natural  variability  of 
the  ecosystem  cannot  be  controlled  (Resh  and  Jackson  1993). 

In  the  same  vein  as  the  step-wise  process  for  identifying  an  ecoregional  biotic 
index  from  raw  metric  values,  the  QC  samples  collected  for  this  project  were  considered 
for  their  precision  and  representativeness  of  the  biological  condition.  With  RBPs  and  the 
use  of  multimetric  assessment  methods,  the  precision  of  the  total  bioassessment  score  is 
as  important  as  the  precision  of  the  individual  metrics  that  comprise  the  score  (Diamond 
et  al.  1996).  Typically,  when  considering  wide  scale  bioassessment  programs,  some  form 
of  criterion  is  established  to  assess  the  quality  of  the  data  that  has  been  collected.  These 
criteria  are  commonly  referred  to  as  Data  Quality  Objectives  (DQOs),  and/or 
Measurement  Quality  Objectives  (MQOs).  These  are  qualitative  and  quantitative 
parameters  developed  by  data  analysts  and  resource  managers  to  evaluate  data  accuracy 
and  quality. 

For  the  Georgia  Ecoregions  Project,  there  was  a  number  of  QA/QC  objectives 
defined  in  the  QAPP.  In  my  research,  only  a  few  of  these  parameters  were  examined. 
Specifically,  analysis  centered  upon  three  measurement  parameters  associated  with  the 
collection  of  benthic  macroinvertebrates:  metric  values,  standardized  metric  scores,  and 
bioassessment  scores.  Not  only  were  these  parameters  analyzed  with  regards  to  the 
prescribed  MQOs,  but  the  QC  data  were  also  evaluated  for  their  relevance  and  value  to 
characterizing  the  reference  and  impaired  conditions,  as  defined  by  the  biotic  index. 

Typical  evaluations  of  data  repeatability  and/or  data  quality  center  on  the 
use  of  a  series  of  calculations  that  quantify  variability  between  measures.  The  following 
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relationships  are  typical  of  bioassessment  protocols  that  define  measures  of  acceptable 
variability  (Stribling  and  Bressler  2004,  Barbour  et  al.  1999,  USEPA  1995): 

•  Relative  Percent  Difference  (RPD):  used  to  quantify  the  proportional  difference 
between  two  measures  as: 

RPD  =  [(Ci  -  C2)/(Ci  +  C2)]  x  100 
Where  Ci  =  the  larger  of  the  two  values  being  compared,  and  C2  =  the  smaller  of 
the  two  values  being  compared  (Berger  et  al.  1996). 

•  Root  Mean  Square  of  Error  (RMSE):  used  as  an  estimate  of  the  standard 
deviation  of  a  group  of  observations.  The  RMSE  is  determined  by  performing  an 
analysis  of  variance  between  duplicate  samples  to  determine  the  mean  square 
error  (MSE)  that  is  representative  of  within  group  variance. 

•  Coefficient  of  Variability  (CV):  calculated  by  expressing  the  standard  deviation 
as  a  percentage  of  the  mean.  The  coefficient  of  variability  for  a  population  was 
calculated  as: 

CV  =  (RMSE/Y)  x  100 
Where  Y  was  the  mean  of  the  dependent  variable  {e.g.  metric  values,  scores,  etc.). 

Values  associated  with  RPDs  and  RMSEs  characterized  some  level  of  precision 
among  the  parameters  being  analyzed.  As  defined  by  the  QAPP  for  the  Georgia 
Ecoregions  Project,  RPDs  for  metric  values,  metric  scores,  and  bioassessment  scores 
were  defined  to  indicate  some  level  of  data  quality  (see  Table  4).  Additionally,  RMSEs 
for  these  same  measurement  parameters  were  developed  here.    For  each  raw  metric 
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value  and  metric  score,  both  RPDs  and  RMSEs  were  calculated,  while  only  the  candidate 
metrics  used  for  the  development  of  the  biotic  index  were  examined.  These  calculations 
provided  some  indication  of  not  only  the  quality  of  the  data  collected,  but  also  acted  as  a 
measure  of  how  representative  the  biological  data  were  to  each  ecoregion.  This 
evaluation  of  the  QC  data  for  this  project  provides  a  framework  for  data  users  and  water 
resource  managers  to  assess  the  reliability  and  inherent  variability  of  the  proposed  biotic 
indices  for  the  state  of  Georgia. 


Table  4.  Precision  Measurement  Quality  Objectives  for  benthic  macroinvertebrates  as 
defined  in  the  Georgia  Ecoregions  Project  Quality  Assurance  Project  Plan  (QAPP)  (see, 
also,  Barbour  et  al.  1999,  USEPA  1995). 


Measurement  Parameter 

Precision  Level 

Metric  Values 

RPD  <  20% 

RMSE  =  TBD* 

Metric  Scores 

RPD  <  5% 

RMSE  =  TBD* 

Bioassessment  Scores 

RPD  <  5% 

RMSE  =  TBD* 

'TBD,  C'to  be  determined");  RMSE  levels  developed  as  a  result  of  this  study. 
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Results 

Precision  calculations  of  the  measurement  parameters  for  benthic 
macroinvertebrate  metrics  are  presented  at  the  primary  ecoregional  level  (refer  to  Table  1 
and  Figure  1  for  ecoregion  and  subecoregion  classifications  and  descriptions).  For  each 
primary  ecoregion,  calculations  of  RPDs  and  RMSEs,  (as  required  by  the  QAPP 
document),  and  CVs  are  provided  as  averages  for  all  QC  samples  collected,  (both 
reference/impaired  and  spatial/temporal),  inclusive  of  their  ecoregional  designation. 
Average  RPDs,  RMSEs,  and  CVs  are  provided  for  each  raw  metric  value,  standardized 
metric  score,  and  the  final  bioassessment  scores  (biotic  indices)  developed  for  each 
ecoregion.  Additionally,  average  RPD,  RMSE,  and  CV  calculations  are  presented 
relative  to  the  specific  metrics  used  in  the  developed  biotic  indices  for  each  primary 
ecoregion  and  subecoregion  (Gore  et  al.  2005).  At  the  subecoregional  level,  average 
RPD,  RMSE,  and  CV  values  for  raw  metric  values,  standardized  metric  scores,  and 
bioassessment  scores  are  provided  in  the  associated  appendices. 

Also  included  in  the  appendices  for  each  ecoregion  are  RPD,  RMSE,  and  CV 
calculations,  (at  the  ecoregion  and  subecoregion  level),  for  each  stream  class  {i.e. 
reference  vs.  impaired),  and  collectively,  for  all  original  and  QC  sites  sampled  within  the 
ecoregional  designation.  Therefore,  the  averages  specific  to  reference  or  impaired 
streams  in  each  ecoregional  designation  only  include  data  relevant  to  the  designated 
stream  class,  whereas  the  averages  provided  for  the  entire  ecoregion  are  inclusive  of  all 
original  and  associated  QC  samples  in  the  ecoregion  designation,  both  reference  and 
impaired,  spatial  and  temporal. 
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Also,  in  conjunction  with  the  required  precision  measures  dictated  by  the  QAPP, 
RPD  values  were  calculated  to  compare  stream  classes,  (i.e.  reference  vs.  impaired),  as 
well  as  QC  sample  type,  (i.e.  spatial  vs.  temporal).  Again,  these  parameters  are 
summarized  for  all  raw  metric  values  at  the  primary  ecoregional  level.  Raw  data 
associated  with  these  parameters,  as  well  as  site  specific  metric  values  for  all  original  and 
QC  sites,  are  provided  on  the  CD-Rom  included  in  the  pocket  materials  of  this  research 
paper. 

RPD  Precision  Measures  for  Raw  Metric  Values  and  Standardized  Metric  Scores 

As  with  determining  the  validity  of  certain  metrics  to  determine  the  biological 
condition  [e.g.  examining  DEs  and  box-and-whisker  plots),  the  raw  metric  values  were 
analyzed  for  RPDs,  acting  both  as  a  measure  of  data  precision  and  data  uncertainty  due  to 
natural  variability  of  the  lotic  ecosystem.  Table  5  contains  a  summary  of  RPDs  averaged 
for  all  metrics  within  each  category.  The  average  RPDs  for  individual  metrics  in  each 
metric  category,  in  most  cases,  were  higher  than  the  measurement  quality  objectives 
dictated  by  the  Ecoregions  QAPP.    The  RPDs  of  raw  metric  values  from  duplicate 


Table  5.    Average  Relative  Percent  Difference  (RPD)  for  all  raw 
each  metric  category  and  per  primary  ecoregion  designation. 

metric  vali 

jes  within 

Metric  Groups 

45 

65 

66 

67 

68 

75 

Taxonomic  Richness 

20.9 

21.4 

20.3 

19.1 

19.2 

18.1 

Community  Composition 

34.2 

32.7 

32.5 

41.5 

27.9 

30.5 

Tolerance/Intolerance 

18.0 

21.3 

19.0 

18.0 

8.0 

18.5 

Functional  Feeding  Group 

25.5 

24.4 

18.5 

26.1 

13.9 

37.2 

Life  Habit 

27.7 

29.2 

25.7 

23.8 

26.9 

28.4 

1 
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reaches  was  expected  to  be  in  80  percent  agreement  (Table  4).  This  is  better  illustrated  in 
Figure  2,  which  demonstrates  that  RPD  values  are  relatively  consistent  between  the 
metric  categories  and  ecoregional  designation. 
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Figure  2.   Average  Relative  Percent  Difference  (RPD)  for  all  raw  metrics  values  within 
each  metric  category  and  per  ecoregion  designation. 


RPDs  were  also  calculated  for  the  standardized  metric  scores  and  have  been 
summarized  for  the  primary  ecoregions  of  Georgia  in  Table  6  and  illustrated  in  Figure  3. 
Again,  these  values  are  higher  than  the  prescribed  MQO  of  95  percent  agreement  for 
standardized  metric  scores. 


Table  6.   Average  Relative  Percent  Difference  (RPD)  for  all  standardized  metric  scores 
within  each  metric  category  and  per  primary  ecoregion  designation. 

Metric  Groups 

45 

65 

66 

67 

68 

75 

Taxonomic  Richness 

18.7 

19.1 

17.8 

16.2 

17.3 

17.3 

Community  Composition 

28.2 

23.8 

28.0 

25.7 

25.2 

20.5 

Tolerance/Intolerance 

23.1 

18.0 

26.1 

11.6 

18.7 

20.7 

Functional  Feeding  Group 

23.4 

22.7 

18.6 

28.1 

8.2 

34.2 

Life  Habit 

27.4 

28.2 

24.0 

24.0 

25.3 

29.4 

Average  RPDs  of  Standarized  Metric  Scores 
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Figure  3.  Average  Relative  Percent  Difference  (RPD)  for  all  standardized  metric  scores 
within  each  metric  category  and  per  primary  ecoregion  designation. 


RAISE  Precision  Measures  of  Raw  Metric  Values  and  Standardized  Metric  Scores 


Another  precision  measurement  utilized  was  the  "root  mean  square  error" 
(RMSE),  which  is  a  representation  of  within  group  variance,  and  acted  as  an  estimate  of 
the  standard  deviation  of  each  population  of  metric  values.  Acceptable  levels  of  error 
associated  with  RMSEs  have  not  been  established  or  quantified  for  this  project.  The 
values  presented  in  this  paper  established  the  ranges  of  variability  on  an  ecoregional  and 
subecoregional  basis.  Similar  to  the  precision  measurements  for  the  RPDs  of  raw  metric 
values  and  standardized  metric  scores,  values  for  RMSEs  presented  here  are  averages  of 
all  metrics  within  each  metric  category  for  each  primary  ecoregion  of  Georgia. 
Subecoregional  averages  of  RMSE  values  are  presented  in  the  associated  appendices. 

The  average  RMSE  was  calculated  for  all  raw  metrics  values  as  summarized  on 
Table  7  and  illustrated  in  Figure  4.    Similarly,  average  RMSEs  for  standardized  metric 
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score  are  presented  in  Table  8  and  illustrated  in  Figure  5. 


Table  7.  Average  Root  Mean  Square  Error  (RMSE)  for  all  raw  metric  values  within  each 
metric  category  and  per  primary  ecoregion  designation. 


Metric  Groups 

45 

65 

66 

67 

68 

75 

Taxonomic  Richness 

1.9 

1.9 

2.7 

1.9 

1.1 

0.9 

Community  Composition 

6.4 

5.9 

5.9 

6.7 

3.7 

5.8 

Tolerance/Intolerance 

3.8 

5.1 

5.1 

5.3 

3.5 

4.4 

Functional  Feeding  Group 

3.5 

3.2 

3.2 

3.3 

2.0 

2.6 

Life  Habit 

2.0 

2.4 

3.2 

2.5 

1.4 

0.9 

Average  RMSEs  of  Raw  Metric  Values 
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Figure  4.  Average  Root  Mean  Square  Error  (RMSE)  for  all  raw  metric  values 
within  each  metric  category  and  per  primary  ecoregion  designation. 


Table  8.  Average  Root  Mean  Square  of  Error  (RMSE)  for  all  standardized  metric  scores 
within  each  metric  category  per  primary  ecoregion  designation. 


Metric  Groups 

45 

65 

66 

67 

68 

75 

Taxonomic  Richness 

10.7 

11.2 

13.2 

13.8 

10.2 

9.4 

Community  Composition 

11.4 

11.2 

13.2 

15.6 

8.8 

10.4 

Tolerance/Intolerance 

10.6 

9.7 

16.2 

11.0 

7.5 

10.6 

Functional  Feeding  Group 

12.3 

12.9 

13.7 

15.3 

6.2 

16.0 

Life  Habit 

12.0 

15.4 

17.6 

15.6 

20.1 

11.5 
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Figure  5.  Average  Root  Mean  Square  Error  (RMSE)  for  all  standardized  metric  scores 
within  each  metric  category  per  primary  ecoregion  designation. 


CV  Precision  Measures  of  Raw  Metric  Values  and  Standardized  Metric  Scores 


The  coefficient  of  variability  (CV)  was  another  measure  of  variability  and 
precision  that  was  calculated  for  raw  metric  values,  standardized  metric  scores,  and 
bioassessment  scores.  Although  not  prescribed  by  the  Georgia  Ecoregions  QAPP  with 
regards  to  MQOs,  CV  values  were  calculated  to  further  illustrate  the  ranges  of  variability 
of  metrics  within  and  between  each  ecoregion.  Statistically,  as  the  CV  value  increases, 
the  precision  of  the  variable  examined  declines.  CV  values  were  calculated  for  raw 
metric  values  (presented  in  Table  9  and  Figure  6),  and  for  standardized  metric  values 
(presented  in  Table  9  and  Figure  7).  Subecoregional  values  of  CVs  for  raw  metric  values 
and  standardized  scores  are  presented  in  the  ecoregional  appendices. 
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Table  9.   Average  Coefficient  of  Variability  (CV)  for  all 
metric  category  per  primary  ecoregion  designation. 

raw  metric  values  within  each 

Metric  Groups 

45 

65 

66 

67 

68 

75 

Taxonomic  Richness 

27.9 

28.0 

28.2 

25.5 

27.1 

28.5 

Community  Composition 

54.5 

55.6 

50.4 

64.4 

39.5 

62.1 

Tolerance/Intolerance 

24.1 

23.3 

26.8 

23.1 

11.3 

26.4 

Functional  Feeding  Group 

32.3 

30.6 

24.2 

27.3 

19.6 

52.8 

Life  Habit 

30.3 

37.6 

37.4 

38.1 

38.1 

41.8 

Average  CVs  of  Raw  Metric  Values 
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Figure  6.  Average  Coefficient  of  Variability  (CV)  for  all  raw  metric  values  within  each 
metric  category  and  per  primary  ecoregion  designation. 


Table  10.  Average  Coefficient  of  Variability  (CV)  for  all  standardized  metric  v 
within  each  metric  category  and  per  primary  ecoregion  designation. 

alues 

Metric  Groups 

45 

65 

66 

67 

68 

75 

Taxonomic  Richness 

23.9 

25.1 

24.0 

24.2 

22.9 

23.7 

Community  Composition 

37.2 

36.3 

33.0 

50.2 

35.7 

29.4 

Tolerance/Intolerance 

25.9 

18.7 

27.0 

15.3 

25.6 

22.6 

Functional  Feeding  Group 

28.1 

28.7 

23.7 

32.4 

11.6 

42.0 

Life  Habit 

28.9 

35.4 

33.9 

27.9 

35.8 

36.2 

35 
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Figure  7.    Average  Coefficient  of  Variability  (CV)  for  all  standardized  metric  scores 
within  each  metric  category  and  per  primary  ecoregion  designation. 


Precision  Measures  for  Bioassessment  Scores 


The  final  precision  measures  for  evaluation  of  the  variability  of  metrics  within 
and  between  ecoregional  designations  centered  upon  the  final  bioassessment  scores  that 
constitute  the  biotic  indices  developed  for  each  ecoregion.  The  average  values  calculated 
for  RPDs,  RMSEs,  and  CVs  presented  in  Table  1 1  were  inclusive  of  only  those  metrics 
that  were  determined  to  be  indicative  of  community  assemblages  that  exhibited  responses 
to  anthropogenic  stress  and  were  descriptive  of  the  reference  and  impaired  condition.  In 
conjunction  with  Table  11,  Figure  8  illustrates  the  ecoregional  averages  of  RPD,  RMSE, 
and  CV  values  for  the  final  bioassessment  scores  used  in  the  development  of  the  biotic 
index.  Comparisons  of  RPD,  RMSE,  and  CV  values  for  final  bioassessment  metrics  for 
each  ecoregion  and  their  corresponding  subecoregions  are  also  presented  in  Tables  12  to 
65  and  illustrated  in  Figures  9  to  13. 
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Table  11.  Average  Relative  Percent  Difference  (RPD),  Root  Mean  Square  Error 
(RMSE),  and  Coefficient  of  Variability  (CV)  values  for  final  bioassessment  scores  used 
in  the  development  of  biotic  indices  for  the  primary  ecoregions  of  Georgia.  (Averages 
are  inclusive  of  only  the  metrics  used  to  develop  the  final  biotic  index  for  each  ecoregion 
designation.) 


RPDs 

45 

65 

66 

67 

68 

75 

20.1 

10.1 

12.7 

10.4 

6.7 

5.4 

RMSEs 

7.6 

5.7 

11.5 

8.6 

6.0 

5.6 

CVs 

22.5 

13.1 

17.2 

14.9 

9.4 

7.7 

Average  Precision  Measurements  for  Final  Bioassessment  Scores 
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Figure  8.  Comparison  of  Average  Relative  Difference  (RPD),  Root  Mean  Square  of 
Error  (RMSE),  and  Coefficient  of  Variability  (CV)  values  for  final  bioassessment  scores 
used  in  the  development  of  biotic  indices  for  the  primary  ecoregions  of  Georgia. 
(Averages  are  inclusive  of  only  the  metrics  used  to  develop  the  final  biotic  index  for  each 
ecoregion  designation.) 

Comparison  of  Precision  Measures  and  Discrimination  Efficiencies  of  Ecoregional 

and  Subecoregional  Biotic  Indices 


The  final  biotic  indices  developed  from  the  Georgia  Ecoregions  study  at  the 
ecoregional  and  subecoregional  level  are  presented  in  the  following  tables.  These  biotic 
indices  were  taken  from  the  Georgia  Ecoregions  numerical  index  report  (Gore  et  al. 
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2005)  submitted  to  Georgia's  DNR  for  establishing  stream  ecosystem  rankings  and 
implementing  guidelines  for  TMDL  permitting  with  the  purpose  of  maintaining  the 
ecological  integrity  of  Georgia's  freshwater  lotic  systems.  In  addition  to  the  following 
tables  that  summarize  the  metrics  used  in  the  biotic  indices  for  the  primary  ecoregions  of 
Georgia,  corresponding  averages  of  the  precision  measures  of  RPDs,  RMSEs,  and  CVs 
are  also  presented  for  the  standardized  metric  scores  that  comprised  the  final  additive 
bioassessment  scores.  Also  included  with  the  precision  measures  are  the  discrimination 
efficiency  values  specific  to  each  metric. 

Summary  of  Biotic  Indices  and  Precision  Measures  for  Ecoregion  45 

Table  12.     The  Biotic  Index  of  primary  Ecoregion  45  developed  from  the  Georgia 
Ecoregions  Project. 

Biotic  Index  -  Ecoregion  45 


Metric 

Metric  Category 

Coleoptera  Taxa 

Richness 

%  Chironomidae 

Composition 

%  Plecoptera 

%  Intolerant  Individuals 

Tolerance  /  Intolerance 

North  Carolina  Biotic  Index  (NCBI) 

Table  13.  Average  precision  measure  values  of  standardized  metric  scores  and 
discrimination  efficiencies  (DE)  for  the  metrics  that  comprise  the  biotic  index  for  the 
primary  ecoregion  45.  (RPD  =  Relative  Percent  Difference;  RMSE  =  Root  Mean  Square 
of  Error;  CV  =  Coefficient  of  Variability) 

Precision  Measures  -  Ecoregion  45 


Metric 

RPD 

RMSE 

CV 

DE 

Coleoptera  Taxa 

32.0 

12.8 

36.1 

0.6 

%  Chironomidae 

34.8 

13.2 

26.7 

0.7 

%  Plecoptera 

29.4 

2.7 

34.0 

0.7 

%  Intolerant  Individuals 

67.0 

18.8 

93.2 

0.6 

North  Carolina  Biotic  Index  (NCBI) 

23.5 

9.9 

17.8 

0.6 

... 
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Table  14.  The  Biotic  Index  of  Subecoregion  45a  developed  from  the  Georgia  Ecoregions 
Project. 

Biotic  Index  -  Subecoregion  45a 


Metric 

Metric  Category 

Plecoptera  Taxa 

Richness 

%  Trichoptera 

Composition 

%  Chironomus  Cricotopus/TC 

Tolerant  Taxa 

Tolerance 

%  Scraper 

Functional  Feeding  Group 

dinger  Taxa 

Habitat 

Table  15.  Average  precision  measure  values  of  standardized  metric  scores  and 
discrimination  efficiencies  (DE)  for  the  metrics  that  comprise  the  biotic  index  for 
Subecoregion  45a.  (RPD  =  Relative  Percent  Difference;  RMSE  =  Root  Mean  Square  of 
Error;  CV  =  Coefficient  of  Variability) 

Precision  Measures  -  Subecoregion  45a 


Metrics 

RPDs 

RMSEs 

CVs 

DEs 

Plecoptera  Taxa 

17.3 

7.6 

37.7 

0.5 

%  Trichoptera 

31.5 

13.4 

30.6 

0.8 

%  Chironomus  Cricotopus/TC 

5.3 

6.4 

6.8 

1.0 

Tolerant  Taxa 

12.2 

10.6 

14.8 

1.0 

%  Scraper 

25.5 

7.6 

22.0 

0.8 

Clinger  Taxa 

10.4 

9.9 

15.4 

0.9 

Table    16.      The  Biotic   Index  of  Subecoregion  45b  developed   from  the   Georgia 
Ecoregions  Project. 

Biotic  Index  -  Subecoregion  45b 


Metric 

Metric  Category 

Coleoptera  Taxa 

Richness 

%  Oligochaeta 

Composition 

%  Plecoptera 

Shredder  Taxa 

Functional  Feeding  Group 

Scraper  Taxa 

Swimmer  Taxa 

Habitat 
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Table  17.  Average  precision  measure  values  of  standardized  metric  scores  and 
discrimination  efficiencies  (DE)  for  the  metrics  that  comprise  the  biotic  index  for 
Subecoregion  45b.  (RPD  =  Relative  Percent  Difference;  RMSE  =  Root  Mean  Square  of 
Error;  CV  =  Coefficient  of  Variability) 

Precision  Measures  -  Subecoregion  45b 


Metrics 

RPDs 

RMSEs 

CVs 

DEs 

Coleoptera  Taxa 

44.4 

8.0 

60.6 

0.9 

%  Oligochaeta 

1.2 

1.4 

1.6 

0.8 

%  Plecoptera 

33.3 

0.6 

141 

0.8 

Shredder  Taxa 

49.2 

10.7 

41.6 

0.9 

Scraper  Taxa 

6.7 

2.7 

10.9 

0.9 

Swimmer  Taxa 

83.3 

31.0 

101 

0.9 

Table  18.  The  Biotic  Index  of  Subecoregion  45c  developed  from  the  Georgia  Ecoregions 
Project. 

Biotic  Index  -  Subecoregion  45c 


Metric 

Metric  Category 

Tanytarsini  Taxa 

Richness 

%  Odonata 

Composition 

%  Tanypodinae/  Total  Chironomidae 

Dominant  Individual 

Tolerance 

%  Intolerant  Individuals 

%  Shredder 

Functional  Feeding  Group 

Swimmer  Taxa 

Habitat 

Table  19.  Average  precision  measure  values  of  standardized  metric  scores  and 
discrimination  efficiencies  (DE)  for  the  metrics  that  comprise  the  biotic  index  for 
Subecoregion  45c.  (RPD  =  Relative  Percent  Difference;  RMSE  =  Root  Mean  Square  of 
Error;  CV  =  Coefficient  of  Variability) 

Precision  Measures  -  Subecoregion  45c 


Metrics 

RPDs 

RMSEs 

CVs 

DEs 

Tanytarsini  Taxa 

0.0 

0.0 

0.0 

0.6 

%  Odonata 

99.8 

35.4 

141 

0.6 

%  Tanypodinae/  TC 

5.4 

6.5 

7.6 

0.6 

Dominant  Individual 

25.6 

8.6 

36.2 

0.6 

%  Intolerant  Individuals 

100 

13.2 

141 

0.8 

%  Shredder 

0.0 

0.0 

0.0 

1.0 

Swimmer  Taxa 

100 

17.7 

141 

0.4 

- 
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Table  20.      The  Biotic  Index  of  Subecoregion  45d  developed  from  the  Georgia 
Ecoregions  Project. 

Biotic  Index  -  Subecoregion  45d 


Metric 

Metric  Category 

Coleoptera  Taxa 

Richness 

%  Tanypodinae/  Total  Chironomidae 

Composition 

%  Odonata 

North  Carolina  Biotic  Index  (NCBI) 

Tolerance 

%  Tolerant  Individuals 

Shredder  Taxa 

Functional  Feeding  Group 

Table  21.  Average  precision  measure  values  of  standardized  metric  scores  and 
discrimination  efficiencies  (DE)  for  the  metrics  that  comprise  the  biotic  index  for 
Subecoregion  45d.  (RPD  =  Relative  Percent  Difference;  RMSE  =  Root  Mean  Square  of 
Error;  CV  =  Coefficient  of  Variability) 

Precision  Measures  -  Subecoregion  45d 


Metrics 

RPDs 

RMSEs 

CVs 

DEs 

Coleoptera  Taxa 

20.0 

8.0 

28.3 

0.8 

%  Tanypodinae/  TC 

0.0 

0.0 

0.0 

1.0 

%  Odonata 

11.8 

14.9 

16.7 

0.8 

NCBI 

5.1 

4.4 

7.2 

1.0 

%  Tolerant  Individuals 

15.9 

17.1 

22.5 

1.0 

Shredder  Taxa 

14.3 

16.5 

20.2 

0.4 

Table   22.      The  Biotic  Index  of  Subecoregion  45h  developed  from  the  Georgia 
Ecoregions  Project. 

Biotic  Index  -  Subecoregion  45h 


Metric 

Metric  Category 

Plecoptera  Taxa 

Richness 

%  Ephemeroptera 

Composition 

%  Plecoptera 

%  Intolerant  Individuals 

Tolerance 

%  Scraper 

Functional  Feeding  Group 

%  dinger 

Habitat 
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Table  23.  Average  precision  measure  values  of  standardized  metric  scores  and 
discrimination  efficiencies  (DE)  for  the  metrics  that  comprise  the  biotic  index  for 
Subecoregion  45h.  (RPD  =  Relative  Percent  Difference;  RMSE  =  Root  Mean  Square  of 
Error;  CV  =  Coefficient  of  Variability) 

Precision  Measures  -  Subecoregion  45h 


Metrics 

RPDs 

RMSEs 

CVs 

DEs 

Plecoptera  Taxa 

18.9 

16.2 

39.7 

1.0 

%  Ephemeroptera 

100 

39.6 

141 

0.8 

%  Plecoptera 

31.0 

20.5 

44.3 

0.8 

%  Intolerant  Individual 

100 

33.4 

141 

0.8 

%  Scraper 

41.0 

25.1 

44.5 

0.6 

%  Clinger 

14.6 

12.4 

21.7 

0.6 

Figure  9  illustrates  the  comparison  of  the  precision  measures  of  RPD,  RMSE,  and 
CV  for  the  biotic  indices  of  Ecoregion  45  and  its  Subecoregions. 
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Figure  9.  Comparison  of  Average  Relative  Difference  (RPD),  Root  Mean  Square  of 
Error  (RMSE),  and  Coefficient  of  Variability  (CV)  values  of  final  metrics  used  in  the 
development  of  biotic  indices  for  Ecoregion  45  and  its  Subecoregions.  (Averages  are 
inclusive  of  only  the  standardized  metrics  used  to  develop  the  final  biotic  index  for  each 
ecoregion  designation.) 
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Summary  of  Biotic  Indices  and  Precision  Measures  for  Ecoregion  65 


Table  24.     The  Biotic  Index  of  primary  Ecoregion  65  developed  from  the  Georgia 
Ecoregions  Project. 

Biotic  Index  -  Ecoregion  65 


Metric 

Metric  Category 

%  Coleoptera 

Composition 

%  Oligochaeta 

Intolerant  Taxa 

Tolerance 

%  Intolerant  Individuals 

%  Predator 

Functional  Feeding  Group 

%  Clinger 

Habitat 

Table  25.  Average  precision  measure  values  of  standardized  metric  scores  and 
discrimination  efficiencies  (DE)  for  the  metrics  that  comprise  the  biotic  index  for  the 
primary  Ecoregion  65.  (RPD  =  Relative  Percent  Difference;  RMSE  =  Root  Mean  Square 
of  Error;  CV  =  Coefficient  of  Variability) 

Precision  Measures  -  Ecoregion  65 


Metric 

RPD 

RMSE 

CV 

DE 

%  Coleoptera 

40.7 

14.1 

43.9 

0.5 

%  Oligochaeta 

5.3 

6.6 

7.1 

0.6 

Intolerant  Taxa 

32.4 

9.3 

23.7 

0.5 

%  Intolerant  Individuals 

42.6 

13.7 

46.2 

0.6 

%  Predator 

23.1 

9.3 

29.8 

0.5 

%  Clinger 

30.3 

14.4 

39.6 

0.5 

Table  26.  The  Biotic  Index  of  Subecoregion  65c  developed  from  the  Georgia  Ecoregions 
Project. 

Biotic  Index  -  Subecoregion  65c 


Metric 

Metric  Category 

%  Trichoptera 

Composition 

Tolerant  Taxa 

Tolerance 

Intolerant  Taxa 

%  Scraper 

Functional  Feeding  Group 

Clinger  Taxa 

Habitat 
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Table  27.  Average  precision  measure  values  of  standardized  metric  scores  and 
discrimination  efficiencies  (DE)  for  the  metrics  that  comprise  the  biotic  index  for 
Subecoregion  65c.  (RPD  =  Relative  Percent  Difference;  RMSE  =  Root  Mean  Square  of 
Error;  CV  =  Coefficient  of  Variability) 

Precision  Measures  -  Subecoregion  65c 


Metrics 

RPDs 

RMSEs 

CVs 

DEs 

%  Trichoptera 

39.5 

25.0 

65.5 

0.7 

Tolerant  Taxa 

6.6 

6.5 

9.0 

0.8 

Intolerant  Taxa 

28.1 

6.7 

12.5 

0.8 

%  Scraper 

58.6 

22.1 

70.6 

0.9 

Clinger  Taxa 

24.3 

19.2 

26.7 

0.6 

Table  28.      The  Biotic  Index  of  Subecoregion  65d  developed  from  the  Georgia 
Ecoregions  Project. 

Biotic  Index  -  Subecoregion  65d 


Metric 

Metric  Category 

Plecoptera  Taxa 

Richness 

%  Chironomidae 

Composition 

%  Hydropsychidae/  EPT 

%  Filterer 

Functional  Feeding  Group 

Swimmer  Taxa 

Habitat 

Table  29.  Average  precision  measure  values  of  standardized  metric  scores  and 
discrimination  efficiencies  (DE)  for  the  metrics  that  comprise  the  biotic  index  for 
Subecoregion  65d.  (RPD  =  Relative  Percent  Difference;  RMSE  =  Root  Mean  Square  of 
Error;  CV  =  Coefficient  of  Variability) 

Precision  Measures  -  Subecoregion  65d 


Metrics 

RPDs 

RMSEs 

CVs 

DEs 

Plecoptera  Taxa 

9.5 

7.9 

12.0 

0.7 

%  Chironomidae 

8.6 

10.0 

12.4 

0.7 

%  Hydropsychidae/  EPT 

13.2 

14.9 

18.6 

0.6 

%  Filterer 

43.1 

30.7 

54.0 

0.7 

Swimmer  Taxa 

12.5 

11.4 

26.0 

0.6 
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Table  30.      The  Biotic  Index  of  Subecoregion  65 g  developed  from  the  Georgia 
Ecoregions  Project. 

Biotic  Index  -  Subecoregion  65g 


Metric 

Metric  Category 

EPT  Taxa 

Richness 

%  Oligochaeta 

Composition 

%  Intolerant  Individuals 

HBI 

Functional  Feeding  Group 

Filterer  Taxa 

Clinger  Taxa 

Habitat 

Table  31.  Average  precision  measure  values  of  standardized  metric  scores  and 
discrimination  efficiencies  (DE)  for  the  metrics  that  comprise  the  biotic  index  for 
Subecoregion  65g.  (RPD  =  Relative  Percent  Difference;  RMSE  =  Root  Mean  Square  of 
Error;  CV  =  Coefficient  of  Variability) 

Precision  Measures  -  Subecoregion  65g 


Metrics 

RPDs 

RMSEs 

CVs 

DEs 

EPT  Taxa 

0.0 

0.0 

0.0 

1.0 

%  Oligochaeta 

100 

18.1 

141 

1.0 

%  Intolerant  Individuals 

0.0 

0.0 

0.0 

1.0 

HBI 

15.2 

4.1 

21.5 

1.0 

Filterer  Taxa 

33.3 

10.6 

47.1 

0.8 

Clinger  Taxa 

20.0 

7.2 

28.3 

1.0 

Table  32.      The  Biotic  Index   of  Subecoregion  65h  developed   from  the  Georgia 
Ecoregions  Project. 

Biotic  Index  -  Subecoregion  65h 


Metric 

Metric  Category 

Tanytarsini  Taxa 

Richness 

Shannon- Wiener  base  e 

%  Oligochaeta 

Composition 

%  Tanytarsini 

NCBI 

Tolerance 

%  Predator 

Functional  Feeding  Group 

Clinger  Taxa 

Habitat 
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Table  33.  Average  precision  measure  values  of  standardized  metric  scores  and 
discrimination  efficiencies  (DE)  for  the  metrics  that  comprise  the  biotic  index  for 
Subecoregion  65h.  (RPD  =  Relative  Percent  Difference;  RMSE  =  Root  Mean  Square  of 
Error;  CV  =  Coefficient  of  Variability) 

Precision  Measures  -  65h 


Metrics 

RPDs 

RMSEs 

CVs 

DEs 

Tanytarsini  Taxa 

64.8 

48.0 

92.7 

0.8 

Shannon- Wiener  base  e 

11.4 

13.4 

15.6 

0.7 

%  Oligochaeta 

7.3 

7.8 

8.9 

0.9 

%  Tanytarsini 

41.4 

10.1 

38.1 

1.0 

/ 

NCBI 

14.5 

14.0 

19.0 

0.8 

%  Predator 

3.4 

2.8 

3.5 

0.6 

dinger  Taxa 

25.1 

23.2 

30.2 

0.9 

Table  34.      The  Biotic   Index  of  Subecoregion  65k  developed  from  the  Georgia 
Ecoregions  Project. 

Biotic  Index  -  Subecoregion  65k 


Metric 

Metric  Category 

%  Gastropoda 

Composition 

%  Orthocladiinae/Total  Chironomidae 

%  Coleoptera 

%  Hydropsychidae/Total  Trichoptera 

%  Filterer 

Functional  Feeding  Group 

%  Collector 

Table  35.      Average  precision  measure  values  of  standardized  metric   scores  and 
discrimination  efficiencies  (DE)  for  the  metrics  that  comprise  the  biotic  index  for 
Subecoregion  65k.   (RPD  =  Relative  Percent  Difference;  RMSE  =  Root  Mean  Square  of 
Error;  CV  =  Coefficient  of  Variability) 
Precision  Measures  -  Subecoregion  65k 


Metrics 


RPDs 


RMSEs 


CVs 


DEs 


%  Gastropoda 
%  Orthocladiinae/Total  Chironomidae 

%  Coleoptera 

%  Hydropsychidae/Total  Trichoptera 

%  Filterer 

%  Collector 


25.0 

14.1 

23.5 

0.8 

63.9 

25.3 

79.3 

0.6 

28.2 

7.2 

24.5 

0.6 

50.0 

35.4 

70.7 

0.6 

27.2 

4.1 

7.5 

0.6 

13.4 

9.5 

16.3 

0.9 

46 


Table  36.  The  Biotic  Index  of  Subecoregion  651  developed  from  the  Georgia  Ecoregions 
Project. 

Biotic  Index  -  Subecoregion  651 


Metric 

Metric  Category 

EPT  Taxa 

Richness 

Diptera  Taxa 

%  EPT 

Composition 

%  Trichoptera 

HBI 

Tolerance 

Predator  Taxa 

Functional  Feeding  Group 

Clinger  Taxa 

Habitat 

Table  37.  Average  precision  measure  values  of  standardized  metric  scores  and 
discrimination  efficiencies  (DE)  for  the  metrics  that  comprise  the  biotic  index  for 
Subecoregion  651.  (RPD  =  Relative  Percent  Difference;  RMSE  =  Root  Mean  Square  of 
Error;  CV  =  Coefficient  of  Variability) 

Precision  Measures  -  651 


Metrics 

RPDs 

RMSEs 

CVs 

DEs 

EPT  Taxa 

75.0 

12.1 

56.6 

0.8 

Diptera  Taxa 

9.4 

9.1 

14.7 

0.6 

%  EPT 

67.0 

3.0 

8.5 

0.8 

%  Trichoptera 

66.7 

7.8 

141 

0.9 

HBI 

19.4 

11.2 

25.5 

0.6 

Predator  Taxa 

8.5 

4.5 

8.3 

0.7 

Clinger  Taxa 

6.7 

4.2 

17.7 

0.8 

Table  38.      The  Biotic   Index  of  Subecoregion  65o   developed   from  the   Georgia 
Ecoregions  Project. 

Biotic  Index  -  Subecoregion  65o 


Metric 

Metric  Category 

Chironomidae  Taxa 

Richness 

%  Oligochaeta 

Composition 

Beck's  Index 

Tolerance 

NCBI 

Scraper  Taxa 

Functional  Feeding  Group 

Burrower  Taxa 

Habitat 

Sprawler  Taxa 
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Table  39.  Average  precision  measure  values  of  standardized  metric  scores  and 
discrimination  efficiencies  (DE)  for  the  metrics  that  comprise  the  biotic  index  for 
Subecoregion  65o.  (RPD  =  Relative  Percent  Difference;  RMSE  =  Root  Mean  Square  of 
Error;  CV  =  Coefficient  of  Variability) 

Precision  Measures  -  Subecoregion  65o 


Metrics 

RPDs 

RMSEs 

CVs 

DEs 

Chironomidae  Taxa 

14.4 

13.5 

19.4 

0.8 

%  Oligochaeta 

3.0 

4.0 

4.2 

0.8 

Beck's  Index 

9.4 

8.1 

10.7 

0.4 

NCBI 

19.2 

16.0 

27.7 

0.6 

Scraper  Taxa 

0.0 

0.0 

0.0 

0.8 

Burrower  Taxa 

18.3 

10.4 

25.9 

0.6 

Sprawler  Taxa 

11.4 

9.8 

11.0 

0.8 

Figure  10  illustrates  the  comparison  of  the  precision  measures  of  RPD,  RMSE, 
and  CV  for  the  biotic  indices  of  Ecoregion  65  and  its  Subecoregions. 
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Figure  10.  Comparison  of  Average  Relative  Difference  (RPD),  Root  Mean  Square  of 
Error  (RMSE),  and  Coefficient  of  Variability  (CV)  values  of  final  metrics  used  in  the 
development  of  biotic  indices  for  Ecoregion  65  and  its  Subecoregions.  (Averages  are 
inclusive  of  only  the  standardized  metrics  used  to  develop  the  final  biotic  index  for  each 
ecoregion  designation. ) 


1 


. 


48 


Summary  of  Biotic  Indices  and  Precision  Measures  for  Ecoregion  66 


Table  40.     The  Biotic  Index  of  primary  Ecoregion  66  developed  from  the  Georgia 
Ecoregions  Project. 

Biotic  Index  -  Ecoregion  66 


Metric 

Metric  Category 

Plecoptera  Taxa 

Richness 

Simpson's  Index 

%  Trichoptera 

Composition 

%  Intolerant  Individuals 

Tolerance 

NCBI 

Predator  Taxa 

Functional  Feeding  Group 

Burrower  Taxa 

Habitat 

Table  41.  Average  precision  measure  values  of  standardized  metric  scores  and 
discrimination  efficiencies  (DE)  for  the  metrics  that  comprise  the  biotic  index  for  primary 
Ecoregion  66.  (RPD  =  Relative  Percent  Difference;  RMSE  =  Root  Mean  Square  of 
Error;  CV  =  Coefficient  of  Variability) 

Precision  Measures  -  Ecoregion  66 


Metric 

RPD 

RMSE 

CV 

DE 

Plecoptera  Taxa 

31.2 

18.8 

33.9 

0.7 

Simpson's  Index 

20.4 

14.0 

20.1 

0.6 

%  Trichoptera 

14.9 

11.5 

16.5 

0.5 

%  Intolerant  Individuals 

24.2 

17.3 

30.1 

0.6 

NCBI 

22.7 

12.8 

20.2 

0.7 

Predator  Taxa 

15.7 

15.7 

21.1 

0.8 

Burrower  Taxa 

21.4 

20.1 

26.6 

0.5 

Table  42.      The  Biotic  Index  of  Subecoregion  66d  developed  from  the  Georgia 
Ecoregions  Project. 

Biotic  Index  -  Subecoregion  66d 


Metric 

Metric  Category 

Diptera  Taxa 

Richness 

%  Plecoptera 

Composition 

%  Odonata 

%  Dominant  Individuals 

Tolerance 

%  Shredder 

Functional  Feeding  Group 

Clinger  Taxa 

Habitat 
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Table  43.  Average  precision  measure  values  of  standardized  metric  scores  and 
discrimination  efficiencies  (DE)  for  the  metrics  that  comprise  the  biotic  index  for 
Subecoregion  66d.  (RPD  =  Relative  Percent  Difference;  RMSE  =  Root  Mean  Square  of 
Error;  CV  =  Coefficient  of  Variability) 

Precision  Measures  -  Subecoregion  66d 


Metrics 

RPDs 

RMSEs 

CVs 

DEs 

Diptera  Taxa 

4.1 

18.5 

26.8 

0.8 

%  Plecoptera 

20.2 

10.8 

24.3 

0.6 

%  Odonata 

50.7 

42.3 

72.6 

1.0 

%  Dominant  Individuals 

3.8 

18.3 

22.5 

0.6 

%  Shredder 

8.6 

15.5 

42.8 

0.8 

dinger  Taxa 

13.8 

16.8 

20.4 

0.6 

Table  44.      The  Biotic  Index  of  Subecoregion  66g  developed  from  the  Georgia 
Ecoregions  Project. 

Biotic  Index  -  Subecoregion  66g 


Metric 

Metric  Category 

EPT  Taxa 

Richness 

%  Chironomidae 

Composition 

%  Tanypodinae/Total  Chironomidae 

NCBI 

Tolerance 

%  Dominant  Individuals 

Scraper  Taxa 

Functional  Feeding  Group 

%  Clinger 

Habitat 

Table  45.  Average  precision  measure  values  of  standardized  metric  scores  and 
discrimination  efficiencies  (DE)  for  the  metrics  that  comprise  the  biotic  index  for 
Subecoregion  66g.  (RPD  =  Relative  Percent  Difference;  RMSE  =  Root  Mean  Square  of 
Error;  CV  =  Coefficient  of  Variability) 

Precision  Measures  -  Subecoregion  66g 


Metrics 

RPDs 

RMSEs 

CVs 

DEs 

EPT  Taxa 

5.4 

10.9 

15.7 

0.9 

%  Chironomidae 

35.9 

15.1 

36.9 

0.9 

%  Tanypodinae/Total  Chironomidae 

35.4 

14.0 

52.3 

0.9 

NCBI 

41.8 

18.4 

28.9 

0.7 

%  Dominant  Individuals 

39.4 

17.0 

25.0 

0.7 

Scraper  Taxa 

13.0 

11.3 

18.8 

0.9 

%  Clinger 

4.9 

10.9 

18.4 

0.7 
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Table  46.  The  Biotic  Index  of  Subecoregion  66j  developed  from  the  Georgia  Ecoregions 
Project. 

Biotic  Index  -  Subecoregion  66j 


Metric 

Metric  Category 

Simpson's  Diversity  Index 

Richness 

Margalef  s  Index 

%  Tanytarsini 

Composition 

%  Intolerant  Individuals 

Tolerance 

Predator  Taxa 

Functional  Feeding  Group 

Sprawler  Taxa 

Habitat 

Table  47.  Average  precision  measure  values  of  standardized  metric  scores  and 
discrimination  efficiencies  (DE)  for  the  metrics  that  comprise  the  biotic  index  for 
Subecoregion  66j.  (RPD  =  Relative  Percent  Difference;  RMSE  =  Root  Mean  Square  of 
Error;  C V  =  Coefficient  of  Variability) 

Precision  Measures  -  Subecoregion  66j 


Metrics 

RPDs 

RMSEs 

CVs 

DEs 

Simpson's  Diversity  Index 

38.9 

26.8 

46.7 

0.8 

Margalef  s  Index 

10.4 

12.1 

14.3 

0.8 

%  Tanytarsini 

63.5 

39.4 

93.0 

0.8 

%  Intolerant  Individuals 

29.9 

22.7 

41.8 

0.6 

Predator  Taxa 

16.1 

14.5 

22.0 

0.8 

Sprawler  Taxa 

19.2 

18.6 

29.7 

0.6 

Figure  1 1  illustrates  the  comparison  of  the  precision  measures  of  RPD,  RMSE, 
and  CV  for  the  biotic  indices  of  Ecoregion  66  and  its  Subecoregions. 
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Figure  11.  Comparison  Average  Relative  Difference  (RPD),  Root  Mean  Square  of  Error 
(RMSE),  and  Coefficient  of  Variability  (CV)  values  of  final  metrics  used  in  the 
development  of  biotic  indices  for  Ecoregion  66  and  its  Subecoregions.  (Averages  are 
inclusive  of  only  the  standardized  metrics  used  to  develop  the  final  biotic  index  for  each 
ecoregion  designation.) 


Summary  of  Biotic  Indices  and  Precision  Measures  for  Ecoregions  67  and  68 


Table  48.     The  Biotic  Index  of  primary  Ecoregion  67  developed  from  the  Georgia 
Ecoregions  Project. 

Biotic  Index  -  Ecoregion  67 


Metric 

Metric  Category 

EPT  Taxa 

Richness 

Plecoptera  Taxa 

%  Plecoptera 

Composition 

%  Isopoda 

HBI 

Tolerance 

Clinger  Taxa 

Habitat 

- 


- 


. 
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Table  49.  Average  precision  measure  values  of  standardized  metric  scores  and 
discrimination  efficiencies  (DE)  for  the  metrics  that  comprise  the  biotic  index  for  primary 
Ecoregion  67.  (RPD  =  Relative  Percent  Difference;  RMSE  =  Root  Mean  Square  of 
Error;  CV  =  Coefficient  of  Variability) 

Precision  Measures  -  Ecoregion  67 


Metric 

RPD 

RMSE 

CV 

DE 

EPT  Taxa 

17.4 

13.5 

24.6 

0.8 

Plecoptera  Taxa 

26.3 

25.7 

61.6 

0.8 

%  Plecoptera 

38.5 

19.9 

130 

0.8 

%  Isopoda 

4.8 

6.4 

7.0 

0.7 

HBI 

8.6 

6.7 

9.6 

0.8 

dinger  Taxa 

7.8 

7.0 

9.9 

0.7 

Table  50.     The  Biotic  Index  of  Subecoregion  67f&i  developed  from  the  Georgia 
Ecoregions  Project. 

Biotic  Index  -  Subecoregion  67f&i 


Metric 

Metric  Category 

EPT  Taxa 

Richness 

Plecoptera  Taxa 

%  EPT 

Composition 

NCBI 

Tolerance 

Scraper  Taxa 

Functional  Feeding  Group 

%  Clinger 

Habitat 

Table  51.  Average  precision  measure  values  of  standardized  metric  scores  and 
discrimination  efficiencies  (DE)  for  the  metrics  that  comprise  the  biotic  index  for 
Subecoregion  67f&i.  (RPD  =  Relative  Percent  Difference;  RMSE  =  Root  Mean  Square 
of  Error;  CV  =  Coefficient  of  Variability) 

Precision  Measures  -  67f&i 


Metrics 

RPDs 

RMSEs 

CVs 

DEs 

EPT  Taxa 

17.1 

16.2 

24.9 

1.0 

Plecoptera  Taxa 

23.1 

25.3 

57.5 

1.0 

%  EPT 

19.7 

17.9 

30.4 

1.0 

NCBI 

7.6 

9.3 

10.3 

0.8 

Scraper  Taxa 

18.3 

18.8 

22.9 

0.8 

%  Clinger 

8.7 

11.7 

15.3 

1.0 

Clinger  Taxa 

10.9 

8.3 

12.1 

0.7 
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Table  52.      The  Biotic  Index  of  Subecoregion  67h  developed  from  the  Georgia 
Ecoregions  Project. 

Biotic  Index  -  Subecoregion  67h 


Metric 

Metric  Category 

Plecoptera  Taxa 

Richness 

%  Gastropoda 

Composition 

%  Tolerant  Individuals 

Tolerance 

HBI 

Scraper  Taxa 

Functional  Feeding  Group 

Swimmer  Taxa 

Habitat 

Table  53.  Average  precision  measure  values  of  standardized  metric  scores  and 
discrimination  efficiencies  (DE)  for  the  metrics  that  comprise  the  biotic  index  for 
Subecoregion  67h.  (RPD  =  Relative  Percent  Difference;  RMSE  =  Root  Mean  Square  of 
Error;  CV  =  Coefficient  of  Variability) 

Precision  Measures  -  67h 


Metrics 

RPDs 

RMSEs 

CVs 

DEs 

Plecoptera  Taxa 

55.6 

47.8 

78.6 

0.5 

%  Gastropoda 

0.0 

0.0 

0.0 

1.0 

%  Tolerant  Individuals 

19.5 

19.7 

27.5 

1.0 

HBI 

0.6 

0.3 

0.9 

1.0 

Scraper  Taxa 

20.0 

12.6 

28.3 

1.0 

Swimmer  Taxa 

33.3 

35.4 

47.1 

1.0 

Table  54.     The  Biotic  Index  of  primary  Ecoregion  68  developed  from  the  Georgia 
Ecoregions  Project. 

Biotic  Index  -  Ecoregion  68 


Metric 

Metric  Category 

Plecoptera  Taxa 

Richness 

%  Hydropsychidae/Total  Trichoptera 

Composition 

%  Tanypodinae/Total  Chironomidae 

NCBI 

Tolerance 

Scraper  Taxa 

Functional  Feeding  Group 

%  Clinger 

Habitat 

• 
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Table  55.  Average  precision  measure  values  of  standardized  metric  scores  and 
discrimination  efficiencies  (DE)  for  the  metrics  that  comprise  the  biotic  index  for  primary 
Ecoregion  68.  (RPD  =  Relative  Percent  Difference;  RMSE  =  Root  Mean  Square  of 
Error;  CV  =  Coefficient  of  Variability) 

Precision  Measures  -  Ecoregion  68 


Metric 

RPD 

RMSE 

CV 

DE 

Plecoptera  Taxa 

0.0 

0.0 

0.0 

0.8 

%  Hydropsychidae/Total  Trichoptera 

7.8 

6.0 

11.1 

0.6 

%  Tanypodinae/Total  Chironomidae 

47.2 

9.0 

66.7 

0.8 

NCBI 

1.8 

2.8 

4.5 

1.0 

Scraper  Taxa 

33.3 

31.7 

47.1 

0.8 

%  Clinger 

1.0 

1.3 

1.4 

0.6 

Figure  12  illustrates  the  comparison  of  the  precision  measures  of  RPD,  RMSE, 
and  CV  for  the  biotic  indices  of  Ecoregion  67  and  its  Subecoregions,  as  well  as 
Ecoregion  68  which  consisted  solely  of  one  Subecoregion  (i.e.  68c&d). 
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Figure  12.  Comparison  of  Average  Relative  Difference  (RPD),  Root  Mean  Square  of 
Error  (RMSE),  and  Coefficient  of  Variability  (CV)  values  of  final  metrics  used  in  the 
development  of  biotic  indices  for  Ecoregion  67  and  its  Subecoregions,  and  Ecoregion  68 
which  consisted  of  one  Subecoregion  (c&d).  (Averages  are  inclusive  of  only  the 
standardized  metrics  used  to  develop  the  final  biotic  index  for  each  ecoregion 
designation.) 
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Summary  of  Biotic  Indices  and  Precision  Measures  for  Ecoregions  75 


Table  56.     The  Biotic  Index  of  primary  Ecoregion  75  developed  from  the  Georgia 
Ecoregions  Project. 

Biotic  Index  -  Ecoregion  75 


Metric 

Metric  Category 

%  Non-Insect 

Composition 

%  Oligochaeta 

%  Odonata 

%  Tanypodinae/Total  Chironomidae 

Tolerant  Taxa 

Tolerance 

%  Tolerant  Individuals 

Table  57.  Average  precision  measure  values  of  standardized  metric  scores  and 
discrimination  efficiencies  (DE)  for  the  metrics  that  comprise  the  biotic  index  for  primary 
Ecoregion  75.  (RPD  =  Relative  Percent  Difference;  RMSE  =  Root  Mean  Square  of 
Error;  CV  =  Coefficient  of  Variability) 

Precision  Measures  -  Ecoregion  75 


Metric 

RPD 

RMSE 

CV 

DE 

%  Non-Insect 

22.2 

9.9 

19.3 

0.6 

%  Oligochaeta 

6.0 

6.2 

6.8 

0.7 

%  Odonata 

2.5 

3.2 

3.4 

0.5 

%  Tanypodinae/Total  Chironomidae 

3.2 

3.4 

3.7 

0.5 

Tolerant  Taxa 

24.1 

17.7 

34.1 

0.6 

%  Tolerant  Individuals 

24.6 

13.8 

26.8 

0.5 

Table  58.  The  Biotic  Index  of  Subecoregion  75e  developed  from  the  Georgia  Ecoregions 
Project. 

Biotic  Index  -  Subecoregion  75e 


Metric 

Metric  Category 

%  Non-Insect 

Composition 

%  Oligochaeta 

%  Isopoda 

%  Odonata 

%  Tolerant  Individuals 

Tolerance 

%  Filterer 

Functional  Feeding  Group 

- 
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Table  59.  Average  precision  measure  values  of  standardized  metric  scores  and 
discrimination  efficiencies  (DE)  for  the  metrics  that  comprise  the  biotic  index  for 
Subecoregion  75e.  (RPD  =  Relative  Percent  Difference;  RMSE  =  Root  Mean  Square  of 
Error;  CV  =  Coefficient  of  Variability) 

Precision  Measures  -  Subecoregion  75e 


Metrics 

RPDs 

RMSEs 

CVs 

DEs 

%  Non- Insect 

18.1 

11.4 

21.4 

0.9 

%  Oligochaeta 

8.2 

7.1 

8.9 

0.9 

%  Isopoda 

28.4 

19.9 

33.3 

0.6 

%  Odonata 

20.2 

21.6 

27.2 

0.6 

%  Tolerant  Individuals 

28.0 

7.2 

28.2 

0.6 

%  Filterer 

41.6 

35.8 

56.1 

0.6 

Table  60.  The  Biotic  Index  of  Subecoregion  75f  developed  from  the  Georgia  Ecoregions 
Project. 

Biotic  Index  -  Subecoregion  75f 


Metric 

Metric  Category 

%  Oligochaeta 

Composition 

%  Tanypodinae/Total  Chironomidae 

Tolerant  Taxa 

Tolerance 

%  Filterer 

Functional  Feeding  Group 

Table  61.  Average  precision  measure  values  of  standardized  metric  scores  and 
discrimination  efficiencies  (DE)  for  the  metrics  that  comprise  the  biotic  index  for 
Subecoregion  75f.  (RPD  =  Relative  Percent  Difference;  RMSE  =  Root  Mean  Square  of 
Error;  CV  =  Coefficient  of  Variability) 

Precision  Measures  -  Subecoregion  75f 


Metrics 

RPDs 

RMSEs 

CVs 

DEs 

%  Oligochaeta 

10.3 

11.6 

14.5 

0.7 

%  Tanypodinae/Total  Chironomidae 

1.9 

2.7 

2.7 

0.8 

Tolerant  Taxa 

50.0 

44.0 

70.7 

0.8 

%  Filterer 

0.0 

0.0 

0.0 

0.7 
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Table  62.      The  Biotic  Index  of  Subecoregion  75h  developed  from  the  Georgia 
Ecoregions  Project. 

Biotic  Index  -  Subecoregion  75h 


Metric 

Metric  Category 

%  Oligochaeta 

Composition 

%  Tolerant  Individuals 

Tolerance 

HBI 

%  Shredder 

Functional  Feeding  Group 

Collector  Taxa 

%  Filterer 

Table  63.  Average  precision  measure  values  of  standardized  metric  scores  and 
discrimination  efficiencies  (DE)  for  the  metrics  that  comprise  the  biotic  index  for 
Subecoregion  75h.  (RPD  =  Relative  Percent  Difference;  RMSE  =  Root  Mean  Square  of 
Error;  CV  =  Coefficient  of  Variability) 

Precision  Measures  -  Subecoregion  75h 


Metrics 

RPDs 

RMSEs 

CVs 

DEs 

%  Oligochaeta 

2.4 

3.4 

3.5 

0.8 

%  Tolerant  Individuals 

9.8 

12.7 

13.9 

0.8 

HBI 

27.7 

29.3 

39.2 

0.5 

%  Shredder 

100 

15.8 

141 

0.8 

Collector  Taxa 

4.3 

4.5 

6.1 

0.5 

%  Filterer 

5.2 

6.5 

7.3 

0.5 

Table  64.  The  Biotic  Index  of  Subecoregion  75j  developed  from  the  Georgia  Ecoregions 
Project. 

Biotic  Index  -  Subecoregion  75j 


Metric 

Metric  Category 

%  Non-Insect 

Composition 

%  Oligochaeta 

%  Tolerant  Individuals 

Tolerance 

Shredder  Taxa 

Functional  Feeding  Group 

• 
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Table  65.  Average  precision  measure  values  of  standardized  metric  scores  and 
discrimination  efficiencies  (DE)  for  the  metrics  that  comprise  the  biotic  index  for 
Subecoregion  75j.  (RPD  =  Relative  Percent  Difference;  RMSE  =  Root  Mean  Square  of 
Error;  CV  =  Coefficient  of  Variability) 

Precision  Measures  -  Subecoregion  75j 


Metrics 

RPDs 

RMSEs 

CVs 

DEs 

%  Non-Insect 

38.0 

37.6 

107 

0.6 

%  Oligochaeta 

7.5 

13.0 

14.2 

0.5 

%  Tolerant  Individuals 

35.4 

34.1 

47.5 

0.6 

Shredder  Taxa 

33.3 

20.4 

49.0 

0.5 

Figure  13  illustrates  the  comparison  of  the  precision  measures  of  RPD,  RMSE, 
and  CV  for  the  biotic  indices  of  Ecoregion  75  and  its  Subecoregions. 
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Figure  13.  Comparison  of  Average  Relative  Difference  (RPD),  Root  Mean  Square  of 
Error  (RMSE),  and  Coefficient  of  Variability  (CV)  values  of  final  metrics  used  in  the 
development  of  biotic  indices  for  Ecoregion  75  and  its  Subecoregions.  (Averages  are 
inclusive  of  only  the  metrics  used  to  develop  the  final  biotic  index  for  each  ecoregion 
designation.) 
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Comparison  of  Relative  Percent  Difference  by  Stream  Class  Designation 

Apart  from  the  requirements  of  the  QAPP  to  analyze  data  quality  and  the 
variability  of  the  biological  metrics  used  in  the  development  of  the  biological  indices  at 
the  ecoregional  and  subecoregional  level,  the  range  of  variability  between  and  among 
stream  classes  was  also  considered.  Although  the  final  determination  of  metrics  that 
represented  the  biological  condition  was  ultimately  based  on  each  metric' s  ability  to 
distinguish  differences  in  the  characteristics  of  a  reference  or  impaired  stream  ecosystem, 
it  was  interesting  to  note  the  measures  of  variability  of  the  stream  classes  themselves. 
Specifically,  the  precision  measure  of  relative  percent  difference  was  considered  to 
illustrate  the  variability  of  raw  metric  values  calculated  for  reference  and  impaired 
streams  separately. 

The  RPD  precision  measures  of  raw  metric  values  for  reference  and  impaired 
stream  classes  are  presented  at  the  ecoregional  level  and  per  metric  group  are 
summarized  in  Table  66.  Additionally,  the  differences  in  variability  between  the  stream 
classes  are  illustrated  in  two  manners:  (1)  each  metric  category  is  compared  individually 
between  ecoregional  designations,  and  (2)  each  metric  category  is  compared  to  each  other 
per  primary  ecoregion  designation.  The  RPD  values  presented  in  Table  66  are  illustrated 
in  Figures  14  to  24.  These  illustrations  demonstrate  variability  for  each  metric  category 
in  relation  to  other  metric  categories,  as  well  as  the  variability  of  each  metric  category 
within  each  ecoregional  designation.  Metric  specific  calculations  of  RPDs  for  each 
stream  class  at  the  ecoregional  and  subecoregional  level  are  included  on  the  CD-Rom  in 
the  pocket  materials  of  this  research  paper. 
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Table  66.  Average  Relative  Percent  Difference  (RPD)  of  Quality  Control  (QC)  sites  per 
stream  class  and  per  primary  ecoregion.  Values  are  averaged  for  all  raw  metric  values 
within  the  metric  group  category,    ("na"  denotes  no  QC  sample  was  collected  for  the 


Metric  Groups 

Class 

45 

65 

66 

67 

68 

75 

Taxonomic  Richness 

Reference 
Impaired 

19.0 
22.6 

22.7 
19.5 

21.3 
18.9 

19.7 
13.0 

na 
19.2 

8.6 
21.9 

Community 
Composition 

Reference 
Impaired 

35.5 
32.0 

31.6 
34.3 

31.6 
33.6 

42.5 

32.5 

na 

27.9 

26.7 
32.0 

Tolerance/Intolerance 

Reference 
Impaired 

16.8 
19.1 

21.0 
21.8 

20.5 
17.0 

18.6 
12.3 

na 
8.0 

20.3 
17.8 

Functional  Feeding 
Group 

Reference 
Impaired 

22.4 
28.1 

26.4 

21.4 

21.4 
14.6 

27.7 
11.9 

na 
13.9 

26.5 
41.4 

Life  Habit 

Reference 
Impaired 

23.8 
31.2 

26.5 
33.1 

29.8 
20.3 

24.7 
15.7 

na 
26.9 

8.6 
36.3 
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Figure  14.  Comparison  of  Relative  Percent  Difference  (RPD)  values  averaged  for  all 
raw  metric  values  of  the  taxonomic  richness  metrics  per  ecoregion  designation  and  for 
stream  class  Quality  Control  (QC)  samples.  (No  Reference  QC  samples  were  collected 
for  Ecoregion  68.) 
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Figure  15.  Comparison  of  Relative  Percent  Difference  (RPD)  values  averaged  for  all 
raw  metric  values  of  the  community  composition  metrics  per  ecoregion  designation  and 
for  stream  class  Quality  Control  (QC)  samples.  (No  Reference  QC  samples  were 
collected  for  Ecoregion  68.) 
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Figure  16.  Comparison  of  Relative  Percent  Difference  (RPD)  values  averaged  for  all 
raw  metric  values  of  the  tolerant/intolerant  individuals  metrics  per  ecoregion  designation 
and  for  stream  class  Quality  Control  (QC)  samples.  (No  Reference  QC  samples  were 
collected  for  Ecoregion  68.) 
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Figure  17.  Comparison  of  Relative  Percent  Difference  (RPD)  values  averaged  for  all 
raw  metric  values  of  the  functional  feeding  group  metrics  per  ecoregion  designation  and 
for  stream  class  Quality  Control  (QC)  samples.  (No  Reference  QC  samples  were 
collected  for  Ecoregion  68.) 
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Figure  18.  Comparison  of  Relative  Percent  Difference  (RPD)  values  averaged  for  all 
raw  metric  values  of  the  life  habit  metrics  per  ecoregion  designation  and  for  stream  class 
Quality  Control  (QC)  samples.  (No  Reference  QC  samples  were  collected  for  Ecoregion 
68.) 
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Figure  19.  Comparison  of  Relative  percent  Difference  (RPD)  values  averaged  for  all 
raw  metric  values  by  stream  class  designation  per  metric  category  for  primary  Ecoregion 
45. 
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Figure  20.  Comparison  of  Relative  percent  Difference  (RPD)  values  averaged  for  all 
raw  metric  values  by  stream  class  designation  per  metric  category  for  primary  Ecoregion 
65. 
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Figure  21. 

raw  metric 
66. 


Comparison  of  Relative  percent  Difference  (RPD)  values  averaged  for  all 
values  by  stream  class  designation  per  metric  category  for  primary  Ecoregion 
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Figure  22.  Comparison  of  Relative  percent  Difference  (RPD)  values  averaged  for  all 
raw  metric  values  by  stream  class  designation  per  metric  category  for  primary  Ecoregion 
67. 
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Figure  23.  Comparison  of  Relative  percent  Difference  (RPD)  values  averaged  for  all 
raw  metric  values  by  stream  class  designation  per  metric  category  for  primary  Ecoregion 
68.  (No  Reference  QC  samples  were  collected  for  Ecoregion  68.) 
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Figure  24.  Comparison  of  Relative  percent  Difference  (RPD)  values  averaged  for  all 
raw  metric  values  by  stream  class  designation  per  metric  category  for  primary  Ecoregion 
75. 
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Comparison  of  Relative  Percent  Difference  for  Spatial  and  Temporal  QC  Samples 

Similar  to  the  analysis  of  average  RPD  of  raw  metric  values  for  all  metrics  within 
each  metric  category  and  per  stream  class  designation,  differences  in  RPDs  between 
spatial  and  temporal  QC  samples  were  also  considered.  RPD  values  are  summarized  in 
Table  67  and  are  also  illustrated  in  corresponding  Figures  25  to  35.  Again,  these  values 
are  illustrated  in  two  manners:  (1)  each  metric  category  is  compared  individually  between 
ecoregional  designations,  and  (2)  each  metric  category  is  compared  to  each  other  per 
primary  ecoregion  designation.  These  illustrations  demonstrate  variability  for  each 
metric  category  in  relation  to  other  metric  categories,  as  well  as  the  variability  of  each 
metric  category  within  each  ecoregional  designation.  Metric  specific  calculations  of 
RPDs  for  each  stream  QC  designation  at  the  ecoregional  and  subecoregional  level  are 
included  on  the  CD-Rom  in  the  pocket  materials  of  this  research  paper. 


Table  67.  Average  Relative  Percent  Difference  (RPD)  of  Quality  Control  (QC)  sites  per 
QC  type  and  per  primary  ecoregion.  Values  are  averaged  for  all  raw  metric  values  within 
the  metric  group  category,  ("na"  denotes  no  QC  sample  was  collected  for  the  QC  type 
designation.) 


Metric  Groups 

QC  Type 

45 

65 

66 

67 

68 

75 

Taxonomic 
Richness 

Spatial 
Temporal 

20.1 
27.4 

17.7 
27.9 

10.3 
21.1 

16.9 
20.6 

19.2 
na 

13.0 

23.5 

Community 
Composition 

Spatial 
Temporal 

32.1 
44.0 

30.5 
36.7 

22.9 
33.0 

30.0 
47.1 

27.9 
na 

19.8 
41.5 

Tolerance/Intolerance 

Spatial 
Temporal 

15.9 
29.9 

21.2 
21.6 

7.4 
20.6 

11.9 

19.4 

8.0 

na 

13.5 
21.3 

Functional  Feeding 
Group 

Spatial 
Temporal 

22.6 
33.7 

22.9 
26.9 

15.8 
19.5 

13.0 
31.4 

13.9 
na 

35.3 
40.4 

Life  Habit 

Spatial 
Temporal 

32.1 
32.7 

36.4 
30.8 

27.2 
35.1 

19.5 
34.5 

32.1 
na 

20.2 
39.7 
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Figure  25.  Comparison  of  Relative  Percent  Difference  (RPD)  values  averaged  for  all 
raw  metric  values  of  the  taxonomic  richness  metrics  per  ecoregion  designation  and  by 
Quality  Control  (QC)  sample  designation.  (No  Temporal  QC  samples  were  collected  for 
Ecoregion  68.) 
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Figure  26.  Comparison  of  Relative  Percent  Difference  (RPD)  values  averaged  for  all 
raw  metric  values  of  the  community  composition  metrics  per  ecoregion  designation  and 
by  Quality  Control  (QC)  sample  designation.  (No  Temporal  QC  samples  were  collected 
for  Ecoregion  68.) 
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Figure  27.  Comparison  of  Relative  Percent  Difference  (RPD)  values  averaged  for  all 
raw  metric  values  of  the  tolerant/intolerant  individuals  metrics  per  ecoregion  designation 
and  by  Quality  Control  (QC)  sample  designation.  (No  Temporal  QC  samples  were 
collected  for  Ecoregion  68.) 
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Figure  28.  Comparison  of  Relative  Percent  Difference  (RPD)  values  averaged  for  all 
raw  metric  values  of  the  functional  feeding  group  metrics  per  ecoregion  designation  and 
by  Quality  Control  (QC)  sample  designation.  (No  Temporal  QC  samples  were  collected 
for  Ecoregion  68.) 
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Figure  29.  Comparison  of  Relative  Percent  Difference  (RPD)  values  averaged  for  all 
raw  metric  values  of  the  life  habit  metrics  per  ecoregion  designation  and  by  Quality 
Control  (QC)  sample  designation.  (No  Temporal  QC  samples  were  collected  for 
Ecoregion  68.) 
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Figure  30.  Comparison  of  Relative  percent  Difference  (RPD)  values  averaged  for  all 
raw  metric  values  by  Quality  Control  (QC)  sample  designation  per  metric  category  for 
primary  Ecoregion  45. 
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Figure  31.  Comparison  of  Relative  percent  Difference  (RPD)  values  averaged  for  all 
raw  metric  values  by  Quality  Control  (QC)  sample  designation  per  metric  category  for 
primary  Ecoregion  65. 
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Figure  32.  Comparison  of  Relative  percent  Difference  (RPD)  values  averaged  for  all 
raw  metric  values  by  Quality  Control  (QC)  sample  designation  per  metric  category  for 
primary  Ecoregion  66. 
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Figure  33.  Comparison  of  Relative  percent  Difference  (RPD)  values  averaged  for  all 
raw  metric  values  by  Quality  Control  (QC)  sample  designation  per  metric  category  for 
primary  Ecoregion  67. 
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Figure  34.  Comparison  of  Relative  percent  Difference  (RPD)  values  averaged  for  all 
raw  metric  values  by  Quality  Control  (QC)  sample  designation  per  metric  category  for 
primary  Ecoregion  68.  (No  Temporal  QC  samples  were  collected  for  Ecoregion  68.) 
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Figure  35.  Comparison  of  Relative  percent  Difference  (RPD)  values  averaged  for  all 
raw  metric  values  by  Quality  Control  (QC)  sample  designation  per  metric  category  for 
primary  Ecoregion  67. 


Discussion  and  Conclusions 

Rapid  bioassessment  is  essentially  a  "biological  shortcut"  in  comparison  with 
impact  assessment  studies  of  the  past,  where  the  goal  is  to  sample  a  wide  range  of  aquatic 
biota  with  the  fastest  methodology  (Metzeling  and  Miller  2001 ).  The  underlying  premise 
of  rapid  bioassessment  is  the  "minimal"  effort  needed  to  characterize  macroinvertebrate 
communities  that  result  in  "maximum"  information.  At  this  point  in  the  evolution  of 
bioassessment  protocols,  minimal  effort  is  reflected  by  the  limited  number  of  sample 
replicates  and  a  limitation  on  the  number  of  collected  organisms  to  be  used  in  metric 
calculations  (Metzeling  and  Miller  2001).  Although  the  stream  conditions  determined  by 
biological  assessments  are  relayed  to  the  general  public  and  water  resource  mangers  as 
narrative  descriptions,  (i.e.  reference  vs.  impaired  lotic  systems,  or  rankings  of 
"good/fair/poor"),  the  final  determination  of  the  biological  condition  is  the  result  of 
quantitative,  numerical  indicators  with  decision  thresholds. 

Measurement  errors  in  an  ecoregional  study,  due  to  its  complexity  and  high  level 
of  effort,  can  be  compounded  from  out-dated  land  use  data,  as  well  as  errors  in  field 
sampling,  laboratory  subsampling,  taxonomic  identification  and  enumeration,  data  entry, 
and  final  metric  calculations.  The  accumulation  of  errors  from  these  multiple  sources 
results  in  uncertainty  and  overall  variability  (Diamond  et  al.  1996;  Clark  and  Whitfield 
1994).  Calculations  of  variance  within  the  biological  parameters  measured  are  necessary 
for  identifying  the  effects  of  measurement  errors  and/or  inherent  differences  between 
sampling  sites  in  relation  to  the  overall  variance  of  a  metric  or  index  on  an  ecoregional 
and  sub-ecoregional  level  (Karr  and  Chu  1997). 
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At  first  glance  it  is  apparent  that  the  majority  of  the  average  RPD  values  for  the 
metric  categories  considered  in  this  study,  (for  both  raw  metric  scores  and  standardized 
metric  scores),  are  above  the  precision  thresholds  of  the  measurement  quality  objectives 
dictated  by  the  QAPP  document  {see  Tables  4  (p.  27),  5  (p.29),  6  (p.30)  and  Figs.  2  and 
3}.  The  raw  metric  values  for  taxonomic  richness  and  tolerance/intolerance  metric 
categories  appear  to  have  better  precision  overall  compared  to  community  composition 
and  life  habit  measures  {Fig.  2  (p.  30)}.  For  each  ecoregion,  the  average  RPD  values  for 
the  taxonomic  richness  and  tolerance/intolerance  metric  categories  fall  closes  to  the 
QAPP  prescribed  MQO  of  20  percent,  while  community  composition  and  life  habit 
measures  are  consistently  above  the  20  percent  precision  threshold. 

After  standardization  of  the  raw  metric  values,  the  average  RPDs  of  the  metric 
scores,  again,  are  above  the  precision  thresholds  established  by  the  QAPP  document  {Fig. 
3  (p.  31)}.  While  the  average  RPDs  of  the  raw  metric  values  for  some  metric  categories 
were  close  to  the  prescribed  precision  threshold  of  twenty  percent,  the  average  RPDs  for 
standardized  metric  scores  for  all  of  the  metric  categories  are  considerably  higher  than 
the  precision  criterion  of  five  percent. 

There  appears  to  be  some  consistency  among  the  averages  of  the  metric  RPDs 
from  one  ecoregion  to  another,  as  well  as  between  the  metric  categories  {Table  5  (p. 
29)}.  With  the  exception  of  RPD  averages  for  metrics  associated  with  measures  of 
tolerance/intolerance  and  taxonomic  richness  which  hovered  at  the  20  percent  precision 
threshold,  the  remaining  RPD  averages  were  typically  much  greater  than  the  prescribed 
MQO.    This  is  better  illustrated  in  Figure  2,  where  it  is  evident  that  metrics  falling  into 
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the  richness  and  tolerance/intolerance  categories  generally  have  less  variability  and  more 
precision  than  the  other  metric  categories. 

As  stated  before,  RMSE  levels  were  to  be  established  as  a  result  of  this  study. 
Ultimately,  the  ranges  of  RMSE  considered  to  be  acceptable  for  a  bioassessment  program 
will  depend  on  the  objectives  of  the  water  resource  manager  and  the  best  professional 
judgment  of  the  data  analyst.  The  overall  measures  of  error  associated  with  biological 
data  determine  some  level  of  data  quality.  Interpretations  of  data  quality  are  important 
for  the  data  user  and  decision  makers  to  evaluate  the  degree  of  the  reliance  on  technical 
and  scientific  information  (Costanza  et  al.  1992).  As  RMSE  values  are  estimates  of  the 
standard  deviation,  which  is  also  considered  as  a  measure  of  precision.  The  assumption 
is  that  the  larger  the  RMSE  value,  the  less  precision,  and/or  greater  variability  within  the 
measures. 

Figure  4  (  p.  32)  illustrated  an  interesting  pattern  of  average  RMSE  values  for  the 
raw  metric  values  associated  with  each  metric  category  in  each  primary  ecoregion.  Each 
metric  category  seemed  to  exhibit  a  "proportional"  trend  in  variability  when  compared  to 
one  another,  with  the  greatest  ranges  of  variability  still  associated  with  the  community 
composition  measures.  After  metric  value  standardization,  the  range  in  variability  for  all 
metric  categories  increased  significantly  when  compared  to  the  RMSE  values  for  the  raw 
metric  values  {see  Figure  5  (p.  33)}.  Additionally,  the  variability  between  each  metric 
category  generally  became  relatively  uniform  within  each  primary  ecoregion. 

For  the  precision  measure  of  CV,  when  compared  to  the  RMSE  values,  an 
opposite  trend  in  variability  ranges  occurred  between  raw  metric  values  and  standardized 
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metric  scores  {see  Figs.  6  (p.  34)  and  7  (p.  35)}.  This  was  most  evident  when  comparing 
variability  ranges  of  the  community  composition  measures  with  CVs  for  the  raw  metric 
values  which  were  consistently  almost  twice  the  value  of  CVs  for  the  standardized  metric 
scores  for  that  metric  category.  CV  values  ranges  for  all  other  metric  categories  did 
decrease  between  raw  metric  values  and  standardized  scores,  but  not  as  dramatically. 
Again,  as  with  most  other  precision  measures  presented,  community  composition 
measures  were  still  consistently  much  higher  in  their  range  of  variability  compared  to  the 
other  metric  categories. 

Table  11  (p.  36)  presented  a  generalized  summary  of  the  precision  measures 
associated  with  the  final  bioassessment  scores  associated  with  each  primary  ecoregion. 
Upon  comparing  the  range  of  variability  of  RPDs,  RMSEs,  and  CVs,  it  should  be  noted 
that  these  average  values  were  inclusive  only  of  the  standardized  metrics  that  were 
included  in  the  final  biotic  index.  Consequently,  these  average  precision  measures  were 
based  on  a  much  smaller  group  of  metrics,  (5  to  8),  than  the  other  values  presented  which 
were  averaged  for  the  entire  suite  of  metrics  within  each  metric  category.  Additionally, 
since  the  metrics  chosen  to  represent  the  biotic  index  for  an  ecoregion  was  based  on  their 
standardized  values,  primarily  those  metric  that  provided  the  highest  additive  score,  and 
exhibiting  the  best  stress  responses,  were  considered  for  the  index.  Although  the 
variability  between  the  metrics  considered  to  quantify  the  final  bioassessment  score  was 
most  likely  minimal,  resulting  in  smaller  ranges  of  variability,  the  MQO  criterion  of  5 
percent  for  RPDs  between  bioassessment  scores  was  still  not  met. 
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When  considering  most  of  the  raw  count  metrics,  (i.e.  number  of  taxa  per  order, 
functional  feeding  groups,  and  habit),  the  values  of  RPDs  are  probably  not  as  important 
as  the  evaluations  of  variance  among  the  final  bioassessment  scores.  Comparing  the  RPD 
values  of  the  raw  metric  values  and  standardized  metric  scores  {Tables  5  and  6  (pp. 
29  and  30)  with  corresponding  Figures  2  and  3},  to  the  RPD  values  for  the  final 
bioassessment  scores  {Table  11  and  Figure  8  (p.  36)},  it  is  apparent  that  the  range  of 
overall  variability  decreased  among  the  measures.  Although  the  RPD  values  for  the  final 
bioassessment  scores,  are  still  higher  than  the  precision  thresholds  dictated  by  the  QAPP 
for  bioassessment  scores,  the  values,  overall,  are  much  closer  to  the  prescribed  five 
percent  threshold,  (as  compared  to  the  RPD  values  for  raw  metric  values  and 
standardized  metric  scores).  As  for  the  other  precision  measures  of  RMSE  and  CV,  the 
range  of  values  associated  with  the  final  bioassessment  scores  are  also  much  narrower 
than  the  RMSEs  and  CVs  calculated  for  raw  metric  values  and  standardized  scores. 

This  initial  examination  of  the  analysis  of  data  precision  leads  to  two  questions: 
(1)  assuming  that  SOP  protocols  for  field  sampling  of  invertebrates  were  followed  with 
minimal  error,  what  are  the  factors  that  could  possibly  influence  the  range  of  variability 
between  an  established  sample  site  and  its  QC  sample?;  and,  (2)  if  we  assume  that  the 
samples  collected  to  determine  the  biological  condition  are  a  valid  or  characteristic 
representation  of  the  ecosystem,  then  how  should  the  data  be  interpreted  in  relationship  to 
the  predetermined  threshold  values  for  precision  and  data  quality? 

Considering  the  sampling  methods  implemented  for  bioassessment  programs, 
there  are  many  factors  to  consider  after  examining  the  RPD  values  for  both  raw  metric 
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values  and  standardized  metric  scores.  The  first  factor  to  consider  is  the  methodology  for 
sampling  the  invertebrate  community.  As  mentioned  before,  some  bioassessment  studies 
have  centered  upon  specific  habitat  types  (i.e.  riffles  and/or  runs)  to  determine  the 
biological  integrity  of  a  lotic  system.  With  the  invertebrate  sampling  protocols  used  for 
this  study,  a  multiple  habitat  sampling  approach,  while  being  more  inclusive  of  the 
assortment  of  macroinvertebrates  in  freshwater  systems,  can  also  lend  itself  to  greater 
variability  because  of  the  mixture  of  habitats  sampled. 

The  "twenty-jab"  method  prescribed  by  the  RBP  protocol  and  the  Georgia  DNR 
(CSU  2000)  was  designed  to  sample  of  a  variety  of  habitats  with  a  relatively  equal  levels 
of  effort  in  proportion  to  those  habitats  that  typically  occur  in  high-  or  low-gradient  lotic 
systems  {see  Table  2  (p.  15)}.  In  instances  where  the  designated  one-  hundred-meter 
sampling  reach  did  not  provide  the  required  distribution  of  effort  among  the  different 
habitat  types,  the  level  of  effort  was  "reallocated"  and  distributed  evenly  among  those 
habitats  that  were  more  dominant  in  the  sampling  reach.  This  may  be  more  of  a  factor 
with  spatial  (200  meter)  QC  samples,  if  the  variety  of  habitats  designated  for  sampling 
did  not  occur  equally  from  one  sampling  reach  to  the  next.  Therefore,  some  raw  metric 
values  may  be  distorted  if  "jabs"  assigned  to  typically  more  productive  habitats  (i.e. 
riffles  and  snags)  are  replaced  by  "jabs"  of  less  productive  habitats  (i.e.  sand),  and  vice- 
versa.  These  changes  in  the  distribution  of  effort  among  habitat  types  can  cause  large 
variations  in  the  invertebrate  assemblages  collected  that  can  ultimately  affect  the  range  of 
metric  values  between  sites  and  within  all  of  the  metric  categories.  In  turn,  this  can 
be  one  factor 
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possibly  responsible  for  higher  and/or  lower  values  RMSEs  and  CVs  than  may  be 
expected,  as  well  as  RPDs  above  the  precision  thresholds  predetermined  for  this  project. 

Another  factor  in  the  inclusion  and/or  exclusion  of  certain  taxa  from  the 
composite  samples  collected  is  the  use  of  random  subsampling.  Caton  (1991)  developed 
a  gridded  screen  technique  to  increase  objectivity  in  the  laboratory  subsampling  of 
benthic  macroinvertebrates.  For  biomonitoring  programs,  subsampling  has  been 
recommended  as  a  valid  and  cost-effective  procedure  where  time  and  monetary  resources 
are  limited.  The  rationale  behind  the  use  of  subsampling  is  twofold,  where  the  level  of 
effort  expended  on  each  sample  collected  is  relatively  equal  and  representative  estimates 
of  the  invertebrate  population  sampled  are  selected  or  picked. 

In  some  instances,  it  is  possible  for  this  sampling  methodology  to  skew  the 
average  values  of  RPDs  and  RMSEs.  In  cases  where  there  may  be  taxa  that  are  rare,  the 
occurrence  of  just  one  organism  picked  from  either  a  primary  or  QC  sample  would  cause 
an  RPD  value  of  "100"  percent  to  be  assigned  to  that  metric  if  there  were  no  occurrence 
of  that  same  taxon  in  the  corresponding  QC  or  primary  sample  from  the  same  stream. 
This  could  be  misleading  when  considering  precision  thresholds  between  samples,  as  the 
samples  are  most  likely  more  similar  than  the  RPD  value  may  indicate.  In  these  instances 
it  may  be  necessary  to  consider  the  raw  taxonomic  data  and  their  effect  on  the  final  value 
of  the  metric  before  assuming  that  the  quality  of  the  data  is  substandard. 

Also,  it  is  important  to  consider  the  characteristics  of  the  metrics  themselves  and 
to  what  ecoregion  and/or  subecoregion  they  are  being  applied.  For  instance,  the  varying 
geomorphology  across  the  state  is  directly  responsible  for  the  variability  of  habitats  and 
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water  chemistry  in  the  lotic  systems  being  analyzed.  This,  in  turn,  will  dictate  the 
presence  and/or  absence  of  certain  taxonomic  groups  and  individual  organisms  based 
upon  habitat  requirements.  Systems  that  are  more  dominated  by  high  gradient,  headwater 
streams  (e.g.  ecoregions  66  and  67)  with  allochthonous  inputs  will  tend  to  have  higher 
percentages  of  shredders  and  scrapers  (e.g.,  Plecoptera  and  Coleoptera),  while 
systems  with  low  gradient  streams  (e.g.  ecoregions  65  and  75)  will  tend  to  have  higher 
numbers  of  filterers  and  collectors  (e.g.  Trichoptera)  (Vannote  et  al.  1980).  Therefore, 
some  metric  calculations  and  corresponding  determinations  of  the  average  RPD,  RMSE, 
and  CV  values  for  an  ecoregion  and/or  subecoregion  may  not  be  truly  indicative  of  the 
range  of  variability  and/or  quality  of  the  data  collected.  Again,  the  data  analyst  may  have 
to  examine  the  individual  sites  within  each  subecoregion  to  asses  the  validity  of  the 
metric  values  for  the  region  that  was  sampled.  Fortunately,  those  candidate  metrics  that 
may  not  be  significant  or  indicative  of  the  invertebrate  assemblages  in  an  ecoregion 
and/or  subecoregion  are  filtered  out  through  DE  calculations  and  box-and-whisker  plot 
examination  during  the  development  of  the  biotic  index. 

In  the  ecoregional  appendices,  subecoregional  and  site-specific  calculations  of 
RPDs  for  raw  metric  values  and  standardized  metric  scores  are  provided  for  each 
subecoregion.  The  initial  interpretation  of  RPDs  for  these  sites  must  be  considered 
cautiously.  In  instances  where  the  RPD  between  two  samples  was  calculated  to  be 
"zero",  there  are  two  scenarios  to  consider:  an  RPD  value  of  "zero"  is  the  result  of  the 
raw  numbers  of  the  metric  for  the  established  sampling  site  and  its  QC  to  be  either  (1) 
equal  in  value,  or  (2)  for  there  to  be  no  occurrence  of  the  organisms  that  define  the  value 
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of  the  metric  in  either  sample.  When  examining  the  MQOs  established  for  this  project,  a 
value  of  "zero"  would  appear  to  indicate  that  the  original  sample  was  a  representative 
sample  and/or  there  was  minimal,  or  no  error  in  performing  the  sampling  (i.e.  high  data 
precision).  However,  the  metrics  themselves  must  be  considered  for  their  ecological 
significance  to  the  target  ecoregion.  Since  they  are  rare  in  that  ecoregion,  an  RPD  of 
"zero"  for  the  number  of  Plecoptera  in  samples  collected  in  the  coastal  plains 
(ecoregion75)  is  not  as  significant  as  RPD  values  for  the  abundant  non-insect  and 
oligochaete  taxa.  The  inclusion  of  certain  metrics  with  minimal  biological  importance  to 
an  ecoregion  or  subecoregion  can  skew  the  overall  average  RPD  values  for  those  regions, 
as  well  as  affect  the  ranges  of  RMSE  and  CV  values. 

For  RPD  values  that  were  calculated  to  be  "100",  indicating  absolute  difference 
between  QC  and  primary  samples,  some  scrutiny  is  deserved,  as  well.  An  RPD  value  of 
100  is  essentially  the  result  of  a  "presence/absence"  scenario,  where  one  sample  may 
have  as  little  as  one  individual,  but  the  corresponding  sample  will  have  no  occurrence  of 
the  same  organism.  Again,  this  value  could  be  misinterpreted,  as  the  presence  of  possibly 
"rare"  individuals  from  one  sample  compared  to  a  corresponding  QC  and/or  primary 
sample  that  may  not  have  the  same  organism  would  not  be  as  significant  as  the  presence 
of  fifty  individuals  where  the  corresponding  sample  may  have  none.  The  data  analyst 
may  need  to  examine  RPD  values  for  individual  sites  to  determine  if  the  presence  or 
absence  of  certain  organisms  is  significant  in  relation  to  the  ecoregion  being  considered. 

After  standardization  of  raw  metric  values  to  produce  metric  scores,  RPDs  and 
RMSEs  are,  again,  calculated  to  determine  precision  estimates  of  the  collected  data. 
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Similarly,  these  values  must  also  be  examined  with  some  scrutiny  as  the  metrics  are  now 
ranked  on  a  similar  scale,  (i.e.  values  are  expressed  on  a  0  to  100  scale),  where  some 
calculated  values  may  be  negative  and  others  may  be  above  the  upper  limit  values  of 
"100".  Those  calculated  metric  scores  which  fell  into  either  one  of  these  categories  were 
changed  to  values  of  "zero"  or  "100"  respectively. 

As  outlined  before,  there  are  many  steps  in  determining  the  appropriate  suite  of 
metrics  that  should  be  used  to  characterized  and  monitor  the  biological  condition  of  an 
ecoregion.  The  metrics  chosen  are  not  only  indicative  of  the  aquatic  assemblages  of  a 
certain  region,  but  are  also  the  most  sensitive  to  anthropogenic  stresses  in  the  ecosystem. 
Tables  1 1  through  22,  (starting  on  p.  36),  display  the  metrics  that  comprised  the  primary 
ecoregion  biotic  indices  and  their  related  precision  measures  for  each  ecoregional 
designation,  as  well  as  their  DE  values.  The  RPD,  RMSE,  and  CV  values  presented  are 
based  upon  the  standardized  metric  scores  for  those  metrics  included  in  the  index. 

With  the  exception  of  a  few  metrics,  the  majority  of  the  RPD  values  for  the 
standardized  metric  scores  of  the  metrics  comprising  each  index  are  relatively  high. 
Considering  the  precision  thresholds  dictated  by  the  QAPP,  the  RPD  values  for 
standardized  metric  scores  should  ideally  be  less  than  or  equal  to  five  percent,  but  in 
relation  to  the  overall  trends  that  RPDs  have  exhibited  for  raw  metric  values  and 
standardized  scores,  these  results  are  consistent.  Conversely,  there  appears  to  be  no 
definitive  correlation  between  the  range  and/or  variability  of  the  precision  measures  and 
the  DE  values  associated  with  each  metric  chosen.  All  DEs  for  the  metrics  that 
formulated  the  final  indices  for  the  primary  ecoregions  did  meet  the  minimal  criteria  of 
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fifty  percent,  but  higher  DE  values  did  not  correspond  with  lower  variability  in  the  RPD, 
RMSE,  or  CV  values  for  those  metrics. 

It  has  been  noted  in  previous  Georgia  Ecoregion  Project  reports  Gore  et  al.  (2004 
and  2005)  that  DE  values  at  the  subecoregional  level  improved,  on  average,  when 
compared  to  the  DE  values  at  the  primary  ecoregional  level.  This  indicated  that  the 
metrics  used  to  characterize  the  biological  condition  at  the  subecoregional  level  were 
more  indicative  of  the  differences  between  the  reference  and  impaired  condition  at  a 
smaller  scale.  Although  DE  values  at  the  subecoregional  level  improved,  there  was  no 
corresponding  trend  in  the  improvement  (i.e.  reduction  in  the  range  of  variability)  of  the 
precision  measures  of  RPD,  RMSE,  and/or  CV  for  those  metrics  that  constitute  the  biotic 
index. 

When  considering  the  ranges  of  variability  between  the  metric  categories  and 
ecoregions  and  examining  the  average  RPD,  RMSE,  and  CV  values  for  the  metric 
categories,  it  must  be  noted  that  the  number  of  QC  sites  among  the  ecoregions  and 
subecoregions  were  not  distributed  evenly.  This  was  primarily  due  to  the  fact  that  QC 
sites  and  those  sites  classified  as  reference  or  impaired  (both  spatially  and  temporally) 
were  randomly  selected.  Therefore,  final  averages  of  the  precision  measurements  for  raw 
metric  values  and  standardized  metric  scores  may  not  have  been  weighted  evenly.  One 
example  of  this  occurs  for  ecoregion  68,  which  is  comprised  of  only  one  subecoregion 
and  had  only  one  QC  sampled  for  the  impaired  stream  class.  Upon  examining  the  overall 
trends  between  the  ecoregions,  and  for  the  majority  of  the  precision  measures  applied, 
ecoregion  68  consistently  exhibited  a  lower  overall  average  RPD,  RMSE,  and  CV  when 
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compared  to  other  ecoregions.  In  instances  where  there  was  a  minimal  number  of  QC 
samples  per  ecoregion  or  subecoregion,  the  RPD,  RMSE,  and  CV  values  associated  with 
data  precision  and  variability  should  be  considered  with  caution,  as  the  number  of 
replicated  samples  may  not  be  sufficient  to  illustrate  ranges  of  variability  within  a  certain 
ecoregion  or  subecoregion. 

Similarly,  the  number  of  QC  sites  designated  as  reference/impaired  and 
spatial/temporal  was  not  distributed  evenly  among  the  stream  class  nor  the  QC  sample 
type.  Although  not  required  by  the  QAPP  for  analysis,  additional  precision  measures 
were  considered  to  examine  possible  influences  of  variability  between  reference  and 
impaired  sites,  as  well  as  spatial  and  temporal  variability  by  metric  category  across  each 
primary  ecoregion.  The  values  presented  in  Table  66  (p.  60)  consisted  of  overall 
averages  of  all  metrics  within  each  metric  group  per  stream  class  and  primary  ecoregional 
designation.  Upon  examination  of  Figures  14  to  24,  {starting  on  p.  60),  there  does  not 
appear  to  be  any  discernable  overall  pattern  between  metric  category  variability  within 
each  ecoregion  (i.e.  there  is  no  consistency  as  to  the  range  of  variability  between 
reference  QCs  when  compared  to  impaired  QCs). 

There  are  only  a  few  exceptions  to  this,  more  specifically  repeatable  differences 
between  stream  classes  in  relation  to  metric  categories  and  ecoregional  designation.  One 
pattern  that  emerged  was  the  difference  in  reference  and  impaired  RPDs  for  the  metric 
categories  of  taxonomic  richness,  functional  feeding  groups,  and  life  habit  from 
ecoregion  75.  These  metric  groups  within  ecoregion  75  have  the  widest  range  of 
variability  between  the  two  stream  classes,  with  the  reference  RPD  values  consistently 
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being  lower.  Additionally,  the  range  of  variability  for  reference  sites  associated  with 
ecoregion  67  are,  for  all  metric  categories,  is  greater  than  the  range  for  the  impaired  sites. 
In  some  respects,  this  particular  pattern  cannot  be  considered  significant  as  the  average 
impaired  values  was  derived  from  only  one  sample. 

Comparisons  of  RPD  values  for  stream  QC  types  were  also  considered  at  the 
primary  ecoregional  level  by  metric  categories.  It  was  apparent  from  the  series  of 
illustrations  {Figures  25  to  35  (starting  on  p.  67)}  derived  from  Table  67  (p.  66),  that 
there  was  consistency  between  the  RPD  values  for  spatial  and  temporal  QC  samples  at 
both  the  ecoregional  and  metric  category  designation.  With  minimal  exceptions, 
temporal  QC  sites  had  higher  ranges  of  variability  when  compared  to  spatial  QC  sites.  In 
general,  this  might  be  an  expected  conclusion  if  the  sites  that  were  originally  chosen  for 
sampling  were  indeed  indicative  of  a  "typical"  reach  of  the  catchment  being  analyzed. 
Also,  considering  that  spatial  QC  samples  were  essentially  collected  at  the  same  time, 
overall  variability  should  be  minimized. 

Looking  at  the  average  RPDs  for  the  temporal  QC  samples,  it  was  evident  that  at 
both  the  ecoregional  level  and  by  metric  category  designation  that  variability  is  much 
greater.  Similarly,  this  might  also  be  an  expected  result  as  there  are  many  factors  that  can 
influence  biological  communities  over  time  (i.e.  rainfall  patterns,  temperature,  etc.).  In 
trying  to  minimize  the  effects  of  temporal  influences,  a  sampling  "index  period",  as 
mentioned  before,  was  utilized  for  all  samples  collected  in  this  project.  For  sites  that 
were  designated  for  temporal  QC  sampling,  field  crews  attempted  to  sample  the  phase 
QC  stream  at  approximately  the  same  time  of  year  as  previously  done. 
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Unfortunately,  a  precise  determination  of  temporal  affects  on  variability  may  not 
be  identifiable  from  one  sampling  period  to  the  next.  One  factor  to  consider  is  the 
number  of  "degree-days"  from  year  to  year  that  cue  the  life  stages  of  freshwater 
macroinvertebrates.  Depending  on  daily  temperature  patterns  between  years,  a  sample 
collected  one  year  may  have  third  or  fourth  instar  nymphs  which  would,  typically,  be 
easier  to  identify  to  a  lower  taxonomic  level  (i.e.  genus  and/or  species  versus  family 
and/or  order).  Other  corresponding  temporal  QC  samples  may  not  have  had  the  same 
number  of  degree-days  before  sampling  that  could  have  resulted  in  earlier  instar  nymphs 
that  may  not  be  identifiable  to  the  same  taxonomic  resolution  as  a  previous  sample. 
Ultimately,  this  will  affect  values  of  metrics  that  require  lower  taxonomic  resolution  to  be 
quantified. 

More  importantly,  another  temporal  factor  to  consider  is  changes  in  land  use 
patterns  across  the  state  of  Georgia.  For  some  ecoregions,  (for  example  ecoregion  45), 
that  have  highly  urbanized  areas,  (i.e.  Atlanta),  land  use  can  change  on  a  weekly  basis. 
Considering  that  most  land  use  data  used  to  identify  reference  and  impaired  catchments 
are  not  updated  on  yearly,  the  variability  of  the  biological  community  from  one  year  to 
the  next  could  be  extreme.  This  can  also  be  exemplified  in  areas  that  have  no  urban 
influence  at  all,  more  specifically,  areas  with  large  amounts  of  acreage  devoted  to 
silviculture.  Although  less  extreme  in  the  rates  of  changes  as  compared  to  urbanized 
areas,  from  one  year  to  the  next  catchments  may  be  either  mature  stands  or  clearcut. 
Even  with  the  required  buffer  strips  emplaced  to  protect  streams  in  these  areas,  there  still 
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exists  the  possibility  of  wide  variances  in  the  biological  community  from  one  sampling 
event  to  the  next. 

Apart  from  the  nuances  of  the  statistical  analysis  of  this  data,  there  are  many 
biological  factors  that  must  also  be  considered  when  interpreting  data  precision  and 
variability.  As  mentioned  previously,  some  consideration  must  be  given  to  individual 
metrics  and  their  applicability  to  the  ecoregion  in  question.  In  some  instances,  the  initial 
analysis  of  raw  metric  values  may  have  not  provided  DE  values  at  fifty  percent  or  greater, 
which  indicated  the  lack  of  differentiation  between  the  biological  character  of  a  reference 
and  impaired  stream.  One  such  region  that  exemplified  this  was  ecoregion  75  (Coastal 
Plains).  The  metrics  that  comprise  the  biotic  index  for  ecoregion  75  {Table  21  (p.  40)} 
do  not  encompass  measures  from  each  metric  category  (i.e.  taxonomic  richness, 
functional  feeding  groups,  etc.),  but  are  limited  to  the  measures  of  community 
composition  and  tolerant  individuals.  This  is  a  direct  reflection  of  the  invertebrate 
assemblages  that  are  most  characteristic  of  the  ecosystem  in  that  region  and  are  dictated, 
in  part,  by  habitat  features,  (i.e.  predominance  of  sand  and  silt  substrates,  presence  of 
woody  debris,  etc.).  The  majority  of  the  organisms  in  the  coastal  plain  ecoregion  is  non- 
insect  (e.g.  amphipoda,  isopoda,  gastropoda,  oligochaeta,  etc.),  being  poorly  described  or 
quantified  by  traditional  richness  metrics  (Gore  et  al.  2004). 

Additionally,  there  were  some  metrics  associated  with  the  "Life  Habit"  category 
that  were  not  considered  in  development  of  the  biotic  indices  for  all  the  ecoregions. 
Specifically,  the  metrics  for  the  life  habit  category  that  were  excluded  for  use  in  the  biotic 
indices  included  percentages  of  burrowers,  climbers,  sprawlers,  and  swimmers  in  the  lotic 


•'■ 


! 


88 
community.  Although  the  EDAS  database  did  provide  calculations  for  these  metrics  in 
question,  there  exists  no  definitive  scientific  literature  to  support  what  type  of  stress 
response  would  be  demonstrated  by  the  organisms  included  in  those  groups,  but  rather 
their  responses  are  inferred  from  other  "lifestyle"  characteristics  (Barbour  et  al.  1999, 
Fore  et  al.  1996,  and  DeShon  1995).  For  example,  those  benthic  macroinvertebrates 
whose  characteristic  life  habit  are  classified  as  sprawlers,  burrowers,  etc.,  are  also 
categorized  by  other  attributes  such  as  a  taxonomic  order  (i.e.  Ephemeroptera,  Plecoptera, 
etc.),  and7or  a  feeding  mechanisms  (i.e.  predators,  shredders,  collectors,  etc.).  Therefore, 
stress  responses  that  have  been  established  at  an  order,  tolerance,  and/or  feeding  level 
have  been  correlated  to  a  similar  stress  response  for  the  organism  in  accord  with  its  life 
habit  characteristics. 

In  conjunction  with  the  physical  character  of  a  habitat,  the  chemical  character  of 
the  lotic  ecosystem  must  be  considered  as  a  possible  source  of  variability  in 
macroinvertebrate  communities.  Many  ecoregions  and  subecoregions  contained  both 
"clearwater"  and  "blackwater"  streams  within.  In  contrast  to  clearwater  streams, 
blackwater  stream  systems  are  typified  by  high  tannin  inputs  (from  terrestrial  organic 
material),  more  acidic  pH  levels,  and  lower  concentrations  of  dissolved  oxygen.  As  a 
result,  the  benthic  communities  that  dominate  these  systems  can  be  significantly  different 
from  clearwater  streams.  The  physical  habitat  of  blackwater  streams  is  characterized  by 
sandy  substrates  and  fine  particulate  organic  matter,  which  serves  as  an  ideal 
environment  for  oligochaetes,  dipteran  taxa,  and  mollusks,  whereas  clearwater  streams  in 
the  same  region  would  be  dominated  more  by  Trichoptera  taxa    and  acid-intolerant 
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Chironomid  taxa  (Meyer  1990).  There  was  no  separation  of  designated  blackwater  and 
clearwater  streams  in  the  analysis  and  it  has  been  suggested  that  clearwater  and 
blackwater  streams  may  need  to  be  categorized  separately  when  developing  biotic  indices 
for  a  specific  region,  as  these  distinctive  invertebrate  communities  may  respond  to 
anthropogenic  stresses  very  differently  (Gore  et  al.  2004). 

Similarly,  there  has  been  discussion  of  establishing  different  suites  of  metrics  for 
streams  in  ecoregion  75  that  empty  into  the  Atlantic  Ocean.  In  the  initial  selection 
process  for  selecting  reference  streams,  there  was  no  determination  of  the  influences  of 
tidal  effects  within  coastal  catchments.  It  became  evident,  after  taxonomic 
identifications,  that  some  samples  contained  invertebrate  communities  indicative  of 
brackish  water  influences.  Although  these  sites  were  not  initially  considered  to  be 
affected  by  estuarine  influxes,  the  presence  of  salt  tolerant,  marine  species,  such  as 
polychaetes  and  crabs,  presented  some  problems  in  defining  the  biotic  indices  for  streams 
influenced  by  tidal  cycles.  The  primary  problem  is  that  the  EDAS  program  does  not 
account  for  metrics  of  marine  species,  which  creates  difficulty  in  developing  biotic 
indices  that  are  truly  characteristic  of  the  integrity  of  brackish  water  systems.  To  remedy 
this  issue,  Gore  et  al.  (2004)  suggested  identification  of  reference  condition  for  both  so- 
called  "inland"  streams  and  "tidal-coastal"  streams  when  there  is  an  occurrence  of  both 
systems  in  an  ecoregion  or  subecoregion. 

One  major  environmental  factor  that  affected  the  majority  of  the  southeastern 
portion  of  the  state  (i.e.,  below  the  "fall  line"),  was  the  occurrence  of  a  sustained  drought 
over  most  of  the  project's  sampling  phases  from  1999  to  2003.    Many  perennial  streams 
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in  this  region  of  the  state  were  estimated  to  be  dry  for  two  years  or  more.  This  situation 
created  two  problems:  (1)  finding  enough  designated  reference  sites  to  satisfy  the  project 
requirements  for  number  of  sites  to  characterize  the  biological  condition,  and  (2)  when 
designated  reference  sites  did  have  water,  in  most  cases,  there  was  no  initial  indication 
that  the  sample  collected  would  be  representative  of  an  "unstressed,"  typical  biological 
community.  Although  Gore  and  Milner  (1990)  have  demonstrated  that  disturbed  lotic 
systems  can  be  recolonized  by  macroinvertebrates  in  as  little  as  fourteen  to  twenty-one 
days,  there  was  no  sure  way  for  field  teams  to  verify  that  normal  stream  functions  and 
invertebrate  communities  had  returned  to  their  typical  character.  In  instances  where  lotic 
systems  have  had  some  form  of  sustained  stress,  additional  sampling  should  be 
performed  so  that  the  reference  and/or  impaired  conditions  can  be  adequately  defined. 

With  all  of  the  possible  influences  of  the  biological  variables  discussed, 
ecological  responses  to  varying  levels  and  types  of  stressors  can  be  complex  and  difficult 
to  accurately  measure  with  a  high  degree  of  reliability  (Murtaugh  1996).  In  some  cases, 
the  use  of  benthic  macroinvertebrates  for  bioassessment  may  not  provide  a  clear  response 
to  anthropogenic  influences.  Floods  and  droughts  inevitably  will  affect  aquatic 
ecosystems  over  the  course  of  time.  Where  there  are  instances  of  "pulse"  events  in  an 
ecosystem,  the  induced  stress  may  not  be  significant  enough  to  permanently  alter  the 
composition  of  the  aquatic  community,  especially  invertebrates  with  their  ability  to 
recolonize  (Gore  and  Milner  1990),  fecundity,  and  dispersal  ability  (Patrick  1975).  Pulse 
events  can  be  either  naturally  occurring  (i.e.  droughts  and  floods),  or  human  induced  (i.e. 
industrial  discharges). 
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In  bioassessment  programs,  it  is  important  to  identify  reference  and  impaired  sites 
that  encompass  natural  variability  within  and  between  watersheds,  as  well  as  the 
variability  of  the  influences  of  possible  anthropogenic  influences.  This  is  crucial  to  the 
subsequent  calculated  metrics  and  biotic  indices  that  water  resource  managers  will  utilize 
in  their  decision-making  process.  Evaluation  of  stream  ecosystem  health  can  be  hindered 
by  the  cost  and  time  constraints  posed  by  large  scale  quantitative  biomonitoring  sampling 
protocols  (Rosenberg  and  Resh  1993).  The  goal  of  quality  control  protocols  is  to 
measure  the  quality  of  a  procedure  so  that  it  meets  the  needs  of  the  user,  while  aiming  to 
produce  data  that  is  dependable,  adequate,  and  economical  (USEPA  1995). 

The  concept  of  Measurements  Quality  Objectives  (MQOs),  (also  referred  to  as 
Data  Quality  Objectives),  in  bioassessment  programs  is  a  useful  tool  in  evaluating  the 
consistency  of  data  and  limiting  variability  and  potential  sources  of  measurement  error 
(Diamond  et  al.  1996).  When  comparing  two  samples  to  determine  a  level  of  precision, 
acceptable  differences  are  typically  predetermined  by  MQOs.  These  requirements  for 
data  quality  should  ideally  be  based  on  prior  knowledge  of  sampling  procedures  and 
measurement  variables  that  are  specific  to  the  region  and/or  ecosystem  being  studied 
(USEPA  1989).  Since  there  were  no  initial,  widespread  biological  characterizations  of 
the  lotic  invertebrate  communities  in  the  state  of  Georgia  prior  to  this  ecoregional  study, 
the  precision  values  of  the  MQOs  stated  for  this  project  {refer  to  Table  4  (p.  27)}  may 
have  been  unrealistic  and  unattainable  for  the  ranges  of  the  natural  biological  condition. 

Considering  the  precision  thresholds  established  for  this  study,  there  may  be  some 
questions  and/or  concerns  over  the  validity,  repeatability,  and  quality  of  the  data  used  to 
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determine  the  biotic  indices  for  the  ecoregions  of  Georgia.  It  is  obvious  that  the  majority 
of  RPDs  of  the  measurements  parameters  of  metric  values  and  standardized  metric  scores 
are  above  an  established  threshold  which  is  presumed  to  be  indicative  of  some  level  of 
acceptable  data  quality.  It  is  also  evident  that  values  for  RMSE  and  CV  are  highly 
variable.  The  initial  interpretation  of  these  results  may  lead  water  resource  managers  to 
believe  that  the  data  are  not  very  precise.  In  reality,  these  precision  thresholds  may  have 
to  be  reevaluated  and  reestablished  by  additional  sampling. 

Because  lotic  invertebrate  communities  can  vary  significantly  between  geographic 
regions,  it  can  be  difficult,  at  first,  to  determine  what  should  typify  the  reference  and/or 
impaired  biological  condition  of  an  ecoregion.  In  some  uses  of  biomonitoring  protocols, 
particularly  the  Index  of  Biotic  Integrity  (IBI)  for  fish  communities,  there  is  some 
expectation  of  what  a  fish  community  should  exhibit  under  a  reference  (or  least 
impacted)  condition  (Karr  et  al.  1986;  USEPA  1990).  It  is  these  reference-expectations 
that  can  hinder  their  application  to  other  geographic  regions  (Simon  and  Lyons,  1995). 
Considering  the  invertebrate  data  produced  by  the  Georgia  Ecoregions  Project,  the 
consistency  of  all  metric  categories  having  average  RPDs  above  the  precision  thresholds 
for  both  raw  metric  values  and  standardized  metric  scores  may  demonstrate  that  the  lotic 
systems  across  the  state  of  Georgia  naturally  have  high  variability  from  year-to-year  and 
spatially  within  the  catchment.  This  in  turn  may  indicate  that  the  established  precision 
thresholds  may  not  be  indicative  of  the  data  quality  for  this  specific  project. 

One  possible  way  to  determine  if  the  variability  within  ecoregions  and 
subecoregions  is  valid,  or  at  least  more  indicative  of  the  ecosystem,  may  be  to  standardize 
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the  number  of  replicate  samples  so  that  there  is  some  equality  in  the  level  of  effort 
expended  for  each  subecoregion.  For  this  study,  a  random  subset  of  sites  was  chosen 
based  on  the  total  number  of  sites  sampled  throughout  Georgia.  As  per  the  requirements 
of  the  QAPP,  ten  percent  of  the  primary  sampling  sites  were  chosen  for  QC  sampling. 
There  was  no  equal  distribution  of  QC  samples  between  or  within  subecoregions  or 
ecoregions.  In  many  scientific  sampling  protocols,  the  premise  behind  the  use  of  an 
equal  level  of  effort  is  to  reduce  bias  and  to  improve  consistency  and  repeatability 
(USEPA  1995;  Plaflan  et  al.  1989). 

If  the  MQOs  for  this  project,  or  any  other  bioassessment  program,  are  not  to  be 
changed  from  the  USEPA  (1995)  guidelines,  then  there  must  be  some  evaluation  of  the 
rapid  bioassessment  protocols  used  and  how  it  may  be  altered  to  achieve  some 
established  criteria  for  data  quality.  In  a  related  research  project,  differences  in 
subsample  sizes  were  analyzed  to  determine  if  the  prescribed  two  hundred  organisms  to 
develop  the  biotic  indices  for  the  state  of  Georgia  were  adequate  to  characterize  the 
necessary  biological  criteria.  Rai  (2005)  found  that  in  most  cases,  (for  the  ecoregions 
analyzed),  that  a  subsample  of  three  hundred  organisms  improved  discrimination 
efficiencies,  and  better  characterized  between  the  reference  and  impaired  condition.  With 
that  preliminary  research  suggesting  that  the  RBP  methods  used  for  Georgia  may  need  to 
be  altered,  then  it  may  be  necessary  to  alter  the  bioassessment  protocol  further  to  achieve 
an  acceptable  level  of  data  precision  as  generally  mandated  by  the  USEPA. 

To  begin  with,  with  all  of  the  biological  and  physical  factors  and  variables 
previously  discussed,  more  QC  samples  may  need  to  be  collected  to  better  illustrate 
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ranges  of  variability  between  raw  metrics  and  final  bioassessment  scores.  Especially  in 
light  of  the  comparisons  between  spatial  and  temporal  QC  samples,  and  in  conjunction 
with  the  possible  influences  of  a  3  year  drought,  temporal  variability  may  need  to  be 
addressed  in  more  detail.  As  mentioned  before,  blind  randomization  of  QC  samples  may 
not  fully  describe  the  ranges  of  variability,  especially  on  a  temporal  scale.  A  more 
systematic  QC  sampling  protocol  should  be  employed  for  both  reference  and  impaired 
stream  classes,  as  well  as  from  year  to  year. 

Another  factor  to  consider  in  relation  to  the  variability  among  the  temporal  QC 
sites  is  climate  patterns,  primarily  degree  days  from  one  ecoregion  to  the  next.  From  a 
geographic  and  geomorphological  context,  ecoregions  68,  67,  and  66,  (because  of  latitude 
and  elevation),  will  typically  have  less  degree  days  than  ecoregion  75  with  a  milder 
climate.  As  mentioned  before,  life  stages  could  vary  according  to  temperatures 
throughout  the  year,  which,  in  turn  could  affect  the  emergence  patterns  of  some 
macroinvertebrates.  In  colder  climes,  one  index  period  may  be  sufficient  to  characterize 
the  biological  community,  while  in  warmer  regions,  multiple  index  periods  may  be 
needed  to  ensure  the  collection  of  mutil-voltine  invertebrate  species. 

The  initial  results  from  the  four  phases  of  this  ecoregion  project  have  been  used  to 
develop  biocriteria  specific  for  the  ecoregions  and  sub-ecoregions  across  the  state  of 
Georgia.  The  biocriteria  used  to  develop  the  biotic  indices  must  be  reviewed  through 
additional  sampling  in  the  future.  Given  the  variability  of  hydrologic  cycles  over  time,  as 
well  as  changing  land  use  patterns  as  a  result  of  urbanization  or  agricultural  practices, 
biocriteria  themselves  will  not  remain  static.    It  is  important  to  identify  spatial  and 
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temporal  variability  in  these  aquatic  systems  so  that  biocriteria  can  be  used  wisely  in 
water  management  decisions. 

The  specific  objective  of  biomonitoring  and  bioassessment  projects  is  to  obtain 
the  information  needed  to  accomplish  the  project  goals  and  uses.  The  ultimate  goal  is  to 
characterize  the  biological  condition  of  lotic  ecosystems  and  determine  which  metrics 
adequately  discriminate  between  levels  of  impairment,  whether  the  impairment  is 
minimal  or  severe.  Biological  metrics  and  biotic  indices  are  used  as  a  gauge  of  the 
biological  condition,  as  well  as  being  indicative  of  some  type  of  response  to 
anthropogenic  stress.  These  measures  of  biological  integrity  are  subject  to  change  over 
the  course  of  a  prescribed  study  and/or  continued  monitoring.  Many  metrics  may 
ultimately  be  revised  and/or  reevaluated  for  their  effectiveness  and  applicability  to  the 
ecoregional  character  and  the  expectations  of  the  water  quality  program  in  question. 
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Appendix 

A.  Land  Use  Characteristics  and  Statistics  of  the  Ecoregions  of  Georp 

;ia. 
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From  Gore  et  al.  2004. 
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Appendix  B.  Water 
Project. 


Chemistry  Parameters  Analyzed  for  the  Georgia  Ecoregions 


Parameter 
Measured 

Type  of  Sample 
Taken 

Method  / 

Instrumentation 

Used 

Range  of  Detection 

Ammonia 
(mg/1  as  N) 

Grab  Sample 

EPA  Method  #350.3 

0.03  to  1400NH3-N/L 

Nitrite 
(mg/1  as  N) 

Grab  Sample 

EPA  Method  #354. 1 

0.01  to  1.0mgNO2- 

N/L 

Nitrate  (as  N) 

Grab  Sample 

EPA  Method  #353.3 

0.01  to  1.0mgNO3- 

N/L 

Total  Phosphorus 
(mg/1  as  P) 

Grab  Sample 

EPA  Method  #365.3 

0.01  tol.2mgP/L 

Copper  (mg/1) 

Grab  Sample 

EPA  Method  #220.1 

low  detection  limit  is 
0.1  ppm 

Iron  (mg/1) 

Grab  Sample 

EPA  Method  #236. 1 

low  detection  limit  is 
0.1  ppm 

Manganese  (mg/1) 

Grab  Sample 

EPA  Method  #243.1 

low  detection  limit  is 
0.1  ppm 

Zinc  (mg/1) 

Grab  Sample 

EPA  Method  #289. 1 

low  detection  limit  is 
0.1  ppm 

Conductivity 
(mS/cm) 

In  situ 
Measurement 

HydroLab  H-20 
probe 

1  to  lOOmS/cm 

Dissolved 
Oxygen (%) 

In  situ 
Measurement 

HydroLab  H-20 
probe 

Oto  100% 

Dissolved 
Oxygen  (mg/1) 

In  situ 
Measurement 

HydroLab  H-20 
probe 

0.2  to  18.8  mg/L 

PH 

In  situ 
Measurement 

HydroLab  H-20 
probe 

0  to  14  units 

Turbidity  (NTU) 

In  situ 
Measurement 

HydroLab  H-20 
probe 

5  to  1000  NTU 

Water 
Temperature  (°C) 

In  situ 
Measurement 

HydroLab  H-20 
probe 

-5  to  50°C 

Alkalinity  (mg/1 
as  CaC03) 

Grab  Sample 

EPA  Method  #310.1 

All  concentration 
ranges  of  alkalinity 

Hardness  (mg/1  as 
CaC03) 

Grab  Sample 

EPA  Method  #130.2 

All  concentration 
ranges  of  hardness 

From  Gore  et  al.  2004. 
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Appendix  C.  Physical  Characterization/Water  Quality  Field  Data  Sheet. 


PHYSICAL  CHARACTERIZATION/WATER  QUALITY  FIELD  DATA  SHEET  (FRONT) 

STREAM  NAME 

LOCATION 

STATION  # 

STREAM  CLASS 

LAT                               LONG 

RIVER  BASIN 

STORET  § 

AGENCY 

INVESTIGATORS 

FORM  COMPLETED  BY 

DATE 

AM      PM 

REASON  FOR  SURVEY 

SITE  LOCATION/MAP 


)raw  a  map  of  the  site  and  indicate  the  areas  sampled 


STREAM 
CHARACTERIZATION 


Subsystem  Classification 

9  Perennial  0  Intermittent 


G  Tidal 


Stream  Type 

6  Coldwater 


0  Warmwater 


WEATHER  CONDITIONS 


Now 

9 

3 

9 

(intermittent) 

9  % 


Past  24  hours 


Has  there  been  a  heavy  rain  in  the  last  7  days? 

storm  (heavy  rain)     0  Yes  0  No 

rain  (steady  rain) 

showers  Air  Temperature °  C 


_%     cloud  cover 
clear/sunny 


Other 


i 


106 


Appendix  C.  Physical  Characterization/Water  Quality  Field  Data  Sheet,  (cont.) 


PHYSICAL  CHARACTERIZATION/WATER  QUALITY  FIELD  DATA  SHEET  (BACK) 


RIPARIAN  ZONE/ 
[NSTREAM  FEATURES 


Predominant  Surrounding  Landuse 

8  Forest  0  Commercial 

8  Field/Pasture  9  Industrial 

9  Agricultural  0  Other 

8  Residential 

Local  Watershed  NPS  Pollution 

9  No  evidence  0  Some  potential  sources 
8  Obvious  sources 


Canopy  Cover 

8  Partly  open  0  Partly  shaded 

High  Water  Mark       m 


Local  Water  Erosion 

0  None        0  Moderate  0  Heavy 

Estimated  Stream  Width 


Estimated  Stream  Depth 

©Riffle _m        0Run 

0Pool 


in 


Velocity 


m/sec 


0  Shaded  Estimated  Reach  Length 

Channelized  0  Yes 

Dam  Present  0  Yes 


0No 
0No 


RIPARIAN  VEGETATION 
(18  meter  buffer) 


Indicate  the  dominant  type  and  record  the  dominant  species  present 

0  Trees  0  Shrubs  0  Grasses  9  Herbaceous 

dominant  species  present _________________________„„ 


AQUATIC  VEGETATION 


Indicate  the  dominant  type  and  record  the  dominant  species  present 

9  Rooted  emergent  0  Rooted  submergent  0  Rooted  floating  0  Free  Floating 

9  Floating  Algae  0  Attached  Algae 

dominant  species  present  


Portion  of  the  reach  with  vegetative  cover 


% 


SEDIMENT/  SUBSTRATE 


Odors 

9  Normal  0  Sewage 

Petroleum 

9  Chemical  0  Anaerobic 

9  Other 


0  None 


Deposits 

0  Sludge        0  Sawdust          0  Paper  fiber 

e 

Sand 

0  Relict  shells               0  Other 

Oils 

9  Absent     8  Slight       9  Moderate      0  Profuse 


Looking  at  stones  which  are  not  deeply  embedded, 
are  the  undersides  black  in  color? 

0  Yes  0  No 


WATER  QUALITY 


Temperature_ 


Specific  Conductance_ 

Dissolved  Oxygen 

PH 

Turbidity 


Water  Odors 

ONormal/None 
0  Petroleum 
0  Fishy 

Water  Surface  Oils 

0  Slick        0  Sheen 
8  None       0  Other 


0  Sewage 
0  Chemical 
0  Other 


0  Globs       0  Flecks 


WQ  Instrument  Used 


Turbidity  (if  not  measured) 

0  Clear       0  Slightly  turbid  0  Turbid 

0  Opaque    0  Water  color  0  Other_ 


INORGANIC  SUBSTRATE  COMPONENTS 

(should  add  up  to  100%) 


ORGANIC  SUBSTRATE  COMPONENTS 

(does  not  necessarily  add  up  to  100%) 


Substrate 
Type 


Diameter 


%  Composition  in 
Sampling  Reach 


Substrate 
Type 


Characteristic 


%  Composition  in  Sampling 
Area 


Bedrock 


Detritus 


sticks,  wood,  coarse  plant 
materials  (CPOM) 


Boulder       >  256  mm  (10") 


Cobble        64-256  mm  (2.5"-10" 


Muck-Mud 


black,  very  fine  orgamc 


Gravel 


2-64mm(0.r-2.5") 


Sand 


0.06-2mm  (gritty) 


Marl 


grey,  shell  fragments 


Silt 


0.004-0.06  mm 


Clav 


<  0.004  mm  (slick) 


From  Columbus  State  University's  Quality  Assurance  Project  Plan  (QAPP  2000). 
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Appendix  D.      Habitat  Assessment  Field  Data  Sheets.   [From  Columbus  State 
University's  Quality  Assurance  Project  Plan  (QAPP  200).] 


HABITAT  ASSESSMENT  FIELD  DATA  SHEET  --  HIGH  GRADDZNT  STREAMS  (FRONT) 

STREAM  NAME 

LOCATION 

STATION  n 

LAT 

LONG 

INVESTIGATORS 

FORM  COMPLETED  BY 

DATE 

REASON  FOR  SURVEY 

Habitat 
Parameter 

Condition 

Category 

Optimal 

Suboptimal 

Marginal 

Poor 

Greater  than  70%  of 

40-70%  mix  of  stable 

20-40%  mix  of  stable 

Less  than  20%  stable 

1.  Epifaunal 

substrate  favorable  for 

habitat;  well-suited  for  full 

habitat:  habitat 

habitat;  lack  of  habitat  is 

Substrate/ 

epifaunal  colonization 

colonization  potential: 

availability  less  than 

obvious;  substrate  unstable 

Available  Cover 

and  fish  cover;  mix  of 

adequate  habitat  for 

desirable;  substrate 

or  lacking. 

snags,  submerged  logs, 

mamtenance  of  populations; 

frequently  disturbed  or 

undercut  banks,  cobble 

presence  of  additional 

removed. 

or  other  stable  habitat 

substrate  in  the  form  of 

and  at  stage  to  allow 

newfall,  but  not  yet  prepared 

full  colonization 

for  colonization  (may  rate  at 

potential  (i.e., 

high  end  of  scale). 

logs/snags  that  are  not 

new  fall  and  not 

SCORE 

transient). 

20      19      18      17       16 

15     14      13     12     11 

10      9       8       7       6 

5     4     3     2      10 

Gravel,  cobble,  and 

Gravel,  cobble,  and  boulder 

Gravel,  cobble,  and 

Gravel,  cobble,  and  boulder 

2.  Embeddedness 

boulder  particles  are  0- 

particles  are  25-50% 

boulder  particles  are  50- 

particles  are  more  than  75% 

25%  surrounded  by 

surrounded  bv  fine 

75%  surrounded  by  fine 

surrounded  by  fine 

fine  sediment. 

sediment. 

sediment. 

sediment. 

Layering  of  cobble 

provides  diversity  of 

SCORE 

niche  space 

20     19     18     17      16 

15     14     13    12     11 

10      9       8       7       6 

5     4     3     2      10 

All  four  velocity/depth 

Only  3  of  the  4  regimes 

Only  2  of  the  4  habitat 

Dominated  by  1  velocity/ 

3.  Velocity/Depth 

regimes  present  (slow- 

present  (if  fast-shallow  is 

regimes  present  (if  fast- 

depth  regime  (usually  slow- 

Regime 

deep,  slow-shallow, 

missing,  score  lower  than  if 

shallow  or  slow-shallow 

deep). 

fast-deep,  fast- 

missing  other  regimes). 

are  missing,  score  low). 

shallow).  (Slow  is  < 

SCORE 

0.3  m/s,  deep  is  >  0.5 

m  ) 

20      19      18      17       16 

15     14      13    12    11 

10      9       8       7       6 

5     4     3     2      10 

Little  or  no 

Some  new  increase  in  bar 

Moderate  deposition  of 

Heavy  deposits  of  fine 

4.  Sediment 

enlargement  of  islands 

formation,  mostly  from 

new  gravel,  sand  or  fine 

material,  increased  bar 

Deposition 

or  point  bars  and  less 

gravel,  sand  or  fine 

sediment  on  old  and  new 

development;  more  than 

than  5%  of  the  bottom 

sediment;  5-30%  of  the 

bars;  30-50%  of  the 

50%  of  the  bottom 

affected  bv  sediment 

bottom  affected;  slight 

bottom  affected; 

changing  frequently;  pools 

deposition. 

deposition  in  pools. 

sediment  deposits  at 
obstructions, 
constrictions,  and  bends; 
moderate  deposition  of 

almost  absent  due  to 
substantial  sediment 
deposition. 

SCORE 

pools  prevalent. 

20     19     18     17      16 

15     14     13    12    11 

10      9       8       7       6 

5     4     3     2      10 

Water  reaches  base  of 

Water  fills  >75%  of  the 

Water  fills  25-75%  of 

Very  little  water  in  channel 

5.  Channel  Flow 

both  lower  banks,  and 

available  channel;  or  <25% 

the  available  channel, 

and  mostly  present  as 

Status 

minimal  amount  of 

of  channel  substrate  is 

and/or  riffle  substrates 

standing  pools. 

channel  substrate  is 

exposed. 

are  mostly  exposed. 

SCORE 

exposed. 

20      19      18      17       16 

15     14     13    12    11 

10      9       8       7       6 

5     4     3     2      10 

f 
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Appendix  D.  Habitat  Assessment  Field  Data  Sheets,  (cont.) 

HABITAT  ASSESSMENT  FIELD  DATA  SHEET  -  HIGH  GRADIENT  STREAMS  (BACK) 


Habitat 
Parameter 

Condition  Category 

Optimal 

Suboptiraal 

Marginal 

Poor 

6.  Channel  Alteration 

Channelization  or 

Some  channelization 

Channelization  may  be 

Banks  shored  with  gabion 

dredging  absent  or 

present,  usually  in 

extensive; 

or  cement;  over  80%  of  the 

minimal;  stream  with 

areas  of  bridge 

embankments  or 

stream  reach  channelized 

normal  pattern. 

abutments,  evidence  of 

shonng  structures 

and  disrupted    Instream 

past  channelization, 

present  on  both  banks; 

habitat  greatly  altered  or 

i.e.,  dredging,  (greater 

and  40  to  80%  of 

removed  entirely. 

than  past  20  yr)  may  be 

stream  reach 

present,  but  recent 

channelized  and 

channelization  is  not 

disrupted. 

SCORE 

present 

20     19     18     17      16 

15     14      13     12     11 

10      9       8       7       6 

5     4     3     2      10 

Occurrence  of  nffies 

Occurrence  of  riffles 

Occasional  riffle  or 

Generally  all  flat  water  or 

7.  Frequency  of  Riffles 

relatively  frequent;  ratio 

infrequent,  distance 

bend;  bottom  contours 

shallow  nffies;  poor 

(or  bends) 

of  distance  between 

between  riffles  divided 

provide  some  habitat; 

habitat;  distance  between 

riffles  divided  by  width 

by  the  width  of  the 

distance  between  nffies 

nffies  divided  by  the  width 

of  the  stream  <7:1 

stream  is  between  7  to 

divided  by  the  width  of 

of  the  stream  is  a  ratio  of 

(generally  5  to  7); 

15. 

the  stream  is  between 

>25 

variety  of  habitat  is  key. 

15  to  25. 

In  streams  where  riffles 

are  continuous, 

placement  of  boulders 

or  other  large,  natural 

obstruction  is 

SCORE 

important 

20      19      18      17       16 

15     14     13     12    11 

10      9       8       7       6 

5     4     3     2      10 

8.  Bank  Stability  (score 

Banks  stable;  evidence 

Moderately  stable; 

Moderately  unstable; 

Unstable;  many  eroded 

each  bank) 

of  erosion  or  bank 

infrequent,  small  areas 

30-60%  of  bank  in 

areas;  "raw"  areas  frequent 

failure  absent  or 

of  erosion  mostly 

reach  has  areas  of 

along  straight  sections  and 

Note:  determine  left  or 

minimal,  little  potential 

healed  over.  5-30%  of 

erosion;  high  erosion 

bends;  obvious  bank 

right  side  by  facing 

for  future  problems. 

bank  in  reach  has  areas 

potential  during  floods. 

sloughing;  60-100%  of 

downstream. 

SCORE        (LB) 
SCORE         (UR1 

<5%  of  bank  affected. 

of  erosion. 

bank  has  erosional  scars. 

Left  Bank      10              V 

8             7             6 

5            4            3 

2            1            0 

RiahtRank   10             9 

8            7            6 

5           4           3 

2            1            0 

More  than  90%  of  the 

70-90%  of  the 

50-70%  of  the 

Less  than  50%  of  the 

9.  Vegetative  Protection 

streambank  surfaces 

streambank  surfaces 

streambank  surfaces 

streambank  surfaces 

(score  each  bank) 

and  immediate  riparian 

covered  by  native 

covered  by  vegetation; 

covered  by  vegetation, 

zone  covered  by  native 

vegetation,  but  one 

disruption  obvious; 

disruption  of  streambank 

vegetation,  including 

class  of  plants  is  not 

patches  of  bare  soil  or 

vegetation  is  very  high; 

trees,  understory 

well-represented; 

closely  cropped 

vegetation  has  been 

shrubs,  or  nonwoody 

disruption  evident  but 

vegetation  common; 

removed  to 

macrophytes; 

not  affecting  full  plant 

less  than  one-half  of 

5  centimeters  or  less  in 

vegetative  disruption 

growth  potential  to  any 

the  potential  plant 

average  stubble  height. 

through  grazmg  or 

great  extent;  more  than 

stubble  height 

mowing  minimal  or  not 

one-half  of  the 

remaining. 

evident;  almost  all 

potential  plant  stubble 

plants  allowed  to  grow 

height  remaining. 

SCORE (LB) 

SCORE        (RB) 

naturally 

Left  Bank      10          9 

8            7            6 

5            4            3 

2            1            0 

Right  Bank   10          9 

8           7           6 

5           4           3 

2            1           0 

Width  of  riparian  zone 

Width  of  riparian  zone 

Width  of  riparian  zone 

Width  of  riparian  zone  <6 

10.  Riparian 

>18  meters;  human 

1 2  - 1 8  meters ;  human 

6-12  meters;  human 

meters:  little  or  no  ripanan 

Vegetative  Zone  Width 

activities  (i.e.,  parking 

activities  have 

activities  have 

vegetation  due  to  human 

(score  each  bank  riparian 

lots,  roadbeds,  clear- 

impacted  zone  only 

impacted  zone  a  great 

activities. 

zone) 

cuts,  lawns,  or  crops) 

minimally. 

deal. 

SCORE (LB) 

SCORE        (RB) 

have  not  impacted  zone. 

Left  Bank     10             9 

8            7            6 

5           4           3 

2            1            0 

Right  Bank  10             9 

8           7           6 

5           4           3 

2            1            0 

Total  Score 


' 


I 


, 


.  i 


'     - 


109 


Appendix  D.  Habitat  Assessment  Field  Data  Sheets,  (cont.) 

HABITAT  ASSESSMENT  FIELD  DATA  SHEET  -  LOW  GRADIENT  STREAMS  (FRONT) 


STREAM  NAME 

LOCATION 

STATION  U 

LAT 

LONG 

INVESTIGATORS 

FORM  COMPLETED  BY 

DATE 

TIME                                  AM     PM 

REASON  FOR  SURVEY 

Habitat 
Parameter 

Condition  Category 

Optimal 

Suboptimal 

Marginal 

Poor 

Greater  than  50%  of 

30-50%  mix  of  stable 

10-30%  mix  of  stable 

Less  than  1 0%  stable 

1.  Epifaunal 

substrate  favorable  for 

habitat;  well-suited  for 

habitat;  habitat 

habitat,  lack  of  habitat  is 

Substrate/ 

epifaunal  colonization 

full  colonization 

availability  less  than 

obvious;  substrate  unstable 

Available  Cover 

and  fish  cover;  mix  of 

potential;  adequate 

desirable;  substrate 

or  lacking. 

snags,  submerged  logs, 

habitat  for  maintenance 

frequently  disturbed  or 

undercut  banks,  cobble 

of  populations;  presence 

removed. 

or  other  stable  habitat 

of  additional  substrate 

and  at  stage  to  allow  full 

m  the  form  of  newfall, 

colonization  potential 

but  not  yet  prepared  for 

(i.e.,  logs/snags  that  are 

colonization  (may  rate 

not  new  fall  and  not 

at  high  end  of  scale). 

SCORE 

transient). 

20      19      18      17       16 

15     14     13     12    11 

10      9       8       7       6 

5     4     3     2     10 

Mixture  of  substrate 

Mixture  of  soft  sand, 

All  mud  or  clay  or  sand 

Hard-pan  clay  or  bedrock; 

2.  Pool  Substrate 

matenals,  with  gravel 

mud,  or  clav;  mud  may 

bottom;  little  or  no  root 

no  root  mat  or  vegetation. 

Characterization 

and  firm  sand  prevalent; 

be  dominant;  some  root 

mat;  no  submerged 

root  mats  and 

mats  and  submerged 

vegetation. 

submerged  vegetation 

vegetation  present. 

SCORE 

common. 

20     19     18     17      16 

15     14     13     12    11 

10      9       8       7       6 

5     4     3     2      10 

Even  mix  of  large- 

Majority  of  pools  large- 

Shallow  pools  much 

Majority  of  pools  small- 

3.  Pool  Variability 

shallow,  large-deep, 
small-shallow,  small- 

deep;  very  few  shallow. 

more  prevalent  than 
deep  pools. 

shallow  or  pools  absent. 

SCORE 

deep  pools  present 

20     19     18     17      16 

15     14     13    12    11 

10      9       8       7       6 

5     4     3     2      10 

Little  or  no  enlargement 

Some  new  increase  in 

Moderate  deposition  of 

Heavy  deposits  of  fine 

4.  Sediment 

of  islands  or  point  bars 

bar  formation,  mostly 

new  gravel,  sand  or  fine 

material,  increased  bar 

Deposition 

and  less  than  <20%  of 

from  gravel,  sand  or 

sediment  on  old  and 

development;  more  than 

the  bottom  affected  bv 

fine  sediment;  20-50% 

new  bars;  50-80%  of  the 

80%  of  the  bottom  changing 

sediment  deposition. 

of  the  bottom  affected; 

bottom  affected; 

frequently;  pools  almost 

slight  deposition  m 

sediment  deposits  at 

absent  due  to  substantial 

pools. 

obstructions, 
constrictions,  and 
bends;  moderate 
deposition  of  pools 

sediment  deposition. 

SCORE 

prevalent. 

20     19     18     17      16 

15     14      13     12     11 

10      9       8       7       6 

5     4     3     2      10 

Water  reaches  base  of 

Water  fills  >75%  of  the 

Water  fills  25-75%  of 

Very  little  water  in  channel 

5.  Channel  Flow 

both  lower  banks,  and 

available  channel;  or 

the  available  channel, 

and  mostly  present  as 

Status 

minimal  amount  of 

<25%  of  channel 

and/or  riffle  substrates 

standing  pools. 

channel  substrate  is 

substrate  is  exposed. 

are  mostly  exposed. 

SCORE 

exposed. 

20     19     18      17      16 

15     14     13     12     11 

10      9       8       7       6 

5     4     3     2      10 

; 
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Appendix  D.  Habitat  Assessment  Field  Data  Sheets,  (cont.) 


HABITAT  ASSESSMENT  FIELD  DATA  SHEET  -  LOW  GRADIENT  STREAMS  (BACK) 

Habitat 
Parameter 

Condition  Category 

Optimal 

Suboptimal 

Marginal 

Poor 

6.  Channel 

Channelization  or 

Some  channelization 

Channelization  may  be 

Banks  shored  with  gabion  or 

Alteration 

dredging  absent  or 

present,  usually  in  areas 

extensive; 

cement;  over  80%  of  the 

minimal;  stream  with 

of  bridge  abutments; 

embankments  or 

stream  reach  channelized 

normal  pattern. 

evidence  of  past 

shoring  structures 

and  disrupted.  Instream 

channelization,  i.e., 

present  on  both  banks; 

habitat  greatly  altered  or 

dredging,  (greater  than 

and  40  to  80%  of 

removed  entirely. 

past  20  yr)  may  be 

stream  reach 

present,  but  recent 

channelized  and 

channelization  is  not 

disrupted. 

SCORE 

present 

20     19     18     17      16 

15     14     13    12    11 

10      9       8       7       6 

5     4     3     2      10 

The  bends  in  the  stream 

The  bends  in  the  stream 

The  bends  in  the  stream 

Channel  straight,  waterway 

7.  Channel 

increase  the  stream 

increase  the  stream 

increase  the  stream 

has  been  channelized  for  a 

Sinuosity 

length  3  to  4  times 

length  2  to  3  times  longer 

length  1  to  2  times 

long  distance. 

longer  than  if  it  was  in 

than  if  it  was  in  a  straight 

longer  than  if  it  was  in  a 

a  straight  line.  (Note  - 

line. 

straight  line. 

channel  braiding  is 

considered  normal  in 

coastal  plains  and  other 

low-lying  areas.  This 

parameter  is  not  easily 

SCORE 

rated  in  these  areas.) 

20     19      18      17      16 

15     14      13     12     11 

10      9       8       7       6 

5     4     3     2      10 

Banks  stable;  evidence 

Moderately  stable; 

Moderately  unstable, 

Unstable;  many  eroded 

8.  Bank  Stability 

of  erosion  or  bank 

infrequent,  small  areas  of 

30-60%  of  bank  in 

areas;  "raw"  areas  frequent 

(score  each  bank) 

failure  absent  or 

erosion  mostly  healed 

reach  has  areas  of 

along  straight  sections  and 

minimal;  little  potential 

over.  5-30%  of  bank  in 

erosion;  high  erosion 

bends;  obvious  bank 

for  future  problems. 

reach  has  areas  of 

potential  during  floods. 

sloughing;  60-100%  of  bank 

SCORE        (LB) 

score       mm 

<5%  of  bank  affected. 

erosion. 

has  erosional  scars. 

Left  Bank            10      9 

8            7            6 

5            4            3 

2             1            0 

RiahtBank          10      9 

8            7            6 

5            4            3 

2            1            0 

More  than  90%  of  the 

70-90%  of  the 

50-70%  of  the 

Less  than  50%  of  the 

9.  Vegetative 

streambank  surfaces 

streambank  surfaces 

streambank  surfaces 

streambank  surfaces  covered 

Protection  (score 

and  immediate  npanan 

covered  by  native 

covered  by  vegetation; 

by  vegetation,  disruption  of 

each  bank) 

zone  covered  by  native 

vegetation,  but  one  class 

disruption  obvious: 

streambank  vegetation  is 

vegetation,  including 

of  plants  is  not  well- 

patches  of  bare  soil  or 

very  high;  vegetation  has 

Note:  determine 

trees,  understory 

represented;  disruption 

closely  cropped 

been  removed  to 

left  or  right  side  by 

shrubs,  or  nonwoody 

evident  but  not  affecting 

vegetation  common; 

5  centimeters  or  less  in 

facing  downstream. 

macrophytes; 

full  plant  growth 

less  than  one-half  of  the 

average  stubble  height. 

vegetative  disruption 

potential  to  any  great 

potential  plant  stubble 

through  grazing  or 

extent;  more  than  one- 

height  remaining. 

mowing  minimal  or  not 

half  of  the  potential  plant 

evident;  almost  all 

stubble  height  remaining. 

plants  allowed  to  grow 

SCORE        (LB) 
SCORE         (RR) 

naturally. 

Left  Bank             10      9 

8            7            6 

5           4           3 

2            1            0 

RiahtBank          10      9 

8            7            6 

5           4           3 

2            1            0 

Width  of  riparian  zone 

Width  of  riparian  zone 

Width  of  riparian  zone 

Width  of  npanan  zone  <6 

10.  Riparian 

>18  meters;  human 

12-18  meters;  human 

6-12  meters;  human 

meters:  little  or  no  npanan 

Vegetative  Zone 

activities  (i.e.,  parking 

activities  have  impacted 

activities  have  impacted 

vegetation  due  to  human 

Width  (score  each 

lots,  roadbeds,  clear- 

zone  only  minimally. 

zone  a  great  deal 

activities. 

bank  riparian  zone) 

cuts,  lawns,  or  crops) 
have  not  impacted 

SCORE        fLB) 
SCORE         CRR) 

zone. 

Left  Rank             10      9 

8            7            6 

5           4           3 

2            1            0 

RiahtRank          10      9 

8            7            6 

5            4            3 

2            1            0 

Total  Score 
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Appendix  E.     Geographic  Maps  of  Quality  Control  Sampling  Sites  of  the  Georgia 
Ecoregions  Project. 

Ecoregion  45  -  Piedmont 
Spatial  Quality  Control  Sampling  Sites 


Figure  36.  Spatial  Quality  Control  Sampling  Sites  for  Ecoregion  45. 


Ecoregion  45  -  Spatial  Quality  Control  Sampling  Site  List 


Ecoregion 

Subecoregion 

Site  ID 

Site  Name 

Impairment 
Status 

45 

45a 

HH-18 
45a  -  59 

Whooping  Creek 
Rottenwood  Creek 

Reference 
Impaired 

45b 

45b  -  203 
45b -217 
45b -291 

South  Fork 

Flint  River 

Proctor  Creek 

Impaired 
Impaired 
Impaired 

45c 

45c -17 

Upton  Creek 

Impaired 

45d 

45d-4 

West  Fork  Pumpkin  vine 
Creek 

Reference-GIS 

45h 

45h-l 

Three  Mile  Creek 

Impaired 
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Appendix  E.     Geographic  Maps  of  Quality  Control  Sampling  Sites  of  the  Georgia 
Ecoregions  Project,  (cont.) 

Ecoregion  45  -  Piedmont 
Temporal  Quality  Control  Sampling  Sites 


HH16 


Figure  37.  Temporal  Quality  Control  Sampling  Sites  for  Ecoregion  45. 


Ecoregion  45  -  Temporal  Quality  Control  Sampling  Site  List 


Ecoregion 

Subecoregion 

Site  ID 

Site  Name 

Impairment 
Status 

HH-  16 

Town  Creek 

Impaired 

45a 

HH-  18 

Whooping  Creek 

Reference 

45 

45a  -  59 

Rottenwood  Creek 

Reference 

45h-2 

Powder  Creek 

Reference 

45h 

45h-9 

Mud  Creek 

Reference 

45h-16 

Williams  Creek 

Reference-GIS 

45h-17 

Barnes  Creek 

Reference 

■' 
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Appendix  E.     Geographic  Maps  of  Quality  Control  Sampling  Sites  of  the  Georgia 
Ecoregions  Project,  {cont.) 


Ecoregion  65  -  Southeastern  Plains 
Spatial  Quality  Control  Sampling  Sites 


Figure  38.  Spatial  Quality  Control  Sampling  Sites  for  Ecoregion  65. 


Ecoregion  65  -  Spatial  Quality  Control  Sampling  Site  List 


Ecoregion     Subecoregion       Site  ID 


Site  Name 


Impairment 
Status 


65 


65c 


65d 


65g 


65k 


65c  -  5  Butler  Creek 

65c  -  80  Lanahassee  Creek 


65d-4 


Sally  Branch 


65g  -  84         Trib.  to  Pachitla  Creek 


65k  -  54 

65k  -  68 

65k- 113 

65k -129 


Maiden  Creek 

Crooked  Creek 

Town  Creek 

Trib.  to 

KinchafooneeCreek 


Impaired 
Reference-GIS 


Reference-GIS 


Impaired 


Reference-GIS 

Reference-GIS 

Impaired 

Impaired 
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Appendix  E.     Geographic  Maps  of  Quality  Control  Sampling  Sites  of  the  Georgia 
Ecoregions  Project,  (cont.) 

Ecoregion  65  -Spatial  Quality  Control  Sampling  Site  List  (cont.) 


E co region     Subecoregion       Site  ID 


Site  Name 


65 


651 


65o 


65o  -  3  Olive  Creek 

65o  -  23  Clyatt  Mill  Creek 

65o  -  25  Tributary  to  New  River 


Impairment 
Status 


651  -  160         Trib.  to  Canochee  River  Impaired 

651  -  379  Red  Bluff  Creek  Reference 

651  -  420  Mill  Branch  Impaired 


Impaired 
Reference-GIS 
Reference-GIS 


Ecoregion  65  -  Southeastern  Plains 
Temporal  Quality  Control  Sampling  Sites 


Figure  39.  Temporal  Quality  Control  Sampling  Sites  for  Ecoregion  65. 


li 


t. 


115 


Appendix  E.     Geographic  Maps  of  Quality  Control  Sampling  Sites  of  the  Georgia 
Ecoregions  Project,  (cont.) 

Ecoregion  65  -  Temporal  Quality  Control  Sampling  Site  List 


Ecoregion     Subecoregion        Site  ID 


Site  Name 


Impairment 
Status 


65c 


65 


65d 


65h 


65o 


HH  -  25             Pine  Knot  Creek               Reference 
65c  -  8 Sweetwater  Creek Impaired 


65d  -  4                Sally  Branch 
65 d  -  14         Hannahatchee  Creek 
65d  -  20 Day  Creek 


65h  -  202 
65h  -  206 


Callahan  Branch 
Shaw  Creek 


65o  -  23 


Reference-GIS 

Reference-GIS 

Impaired 


Reference-GIS 
Reference-GIS 


Clyatt  Mill  Creek  Reference-GIS 


Ecoregion  66  -  Blue  Ridge 
Spatial  Quality  Control  Sampling  Sites 


Figure  40.  Spatial  Quality  Control  Sampling  Sites  for  Ecoregion  66. 
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Appendix  E.     Geographic  Maps  of  Quality  Control  Sampling  Sites  of  the  Georgia 
Ecoregions  Project,  (cont.) 


Eeoregion  66  -  Spatial  Quality  Control  Sampling  Site  List 


Ecoreeion 

Subeco  region 

Site  ID 

Site  Name 

Impairment 
Status 

66 

66d 

66d  -  38 

West  Fork  Wolf  Creek 

Impaired 

66g 

66g  -  42 

Trib.  to  Talking  Rock  Creek 

Impaired 

Eeoregion  66  -  Blue  Ridge 
Temporal  Quality  Control  Sampling  Sites 


Figure  41.  Temporal  Quality  Control  Sampling  Sites  for  Eeoregion  66. 


■ 


' 


117 


Appendix  E.     Geographic  Maps  of  Quality  Control  Sampling  Sites  of  the  Georgia 
Ecoregions  Project,  (cont.) 

Ecoregion  66  -  Temporal  Quality  Control  Sampling  Site  List 


Ecoregion 


66 


Subecoregion 

66d 


66g 


66j 


Site  ID  Site  Name 

66d  -  40         Chattahoochee  River 


66g-6 
66g-71 


Holly  Creek 
Yellow  Creek 


66j  -  19             Hothouse  Creek 
66j  -  27           Young  Cane  Creek 
66j  -  2 1 1 Bryan  Creek 


Impairment 
Status 

Reference-GIS 


Reference-GIS 
Impaired 


Reference 

Impaired 

Reference-GIS 


Ecoregion  67  -  Ridge  and  Valley 
Spatial  Quality  Control  Sampling  Sites 


Figure  42.  Spatial  Quality  Control  Sampling  Sites  for  Ecoregion  67. 
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Appendix  E.     Geographic  Maps  of  Quality  Control  Sampling  Sites  of  the  Georgia 
Ecoregions  Project,  (cont.) 


Ecoregion  67  - 

-  Spatial  Quality 

Control  Sampling  Site  List 

Ecoregion    Subecoregion 

Site  ID 

Site  Name                ImmmSSl 

Status 

67                 6?f&1 

67f&i - 16 
67f&i-17 

Cane  Creek               Reference-GIS 
Armuchee  Creek          Reference-GIS 

67h 

67h-5 

Trib.  to  Ruff  Creek              Impaired 

Ecoregion  67  -  Ridge  and  Valley 
Temporal  Quality  Control  Sampling  Sites 


Figure  43.  Temporal  Quality  Control  Sampling  Sites  for  Ecoregion  67. 
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Appendix  E.     Geographic  Maps  of  Quality  Control  Sampling  Sites  of  the  Georgia 
Ecoregions  Project,  (cont.) 


Ecoregion  67  -  Temporal  Quality  Control  Sampling  Site  List 


Ecoregion     Subecoregion  Site  ID 


Site  Name 


Impairment 
Status 


67 


67f&i 


67f&i -  16 
67f&i -  25 
67f&i -  27 


Cane  Creek 
Clarks  Creek 
Dykes  Creek 


Reference-GIS 
Reference-GIS 
Reference-GIS 


Ecoregion  68  -  Southwestern  Appalachians 
Spatial  Quality  Control  Sampling  Sites 


Figure  44.  Spatial  Quality  Control  Sampling  Sites  for  Ecoregion  68. 


Ecoregion  68  -  Spatial  Quality  Control  Sampling  Site  List 


Ecoregion     Subecoregion         Site  ID 


Site  Name 


Stream  Class 


68 


68c&d  68c&d-8        Tributary  to  Middle  Fork         Impaired 

Little  River 


! 

: 


■ 


. 


■ 


120 

Appendix  E.     Geographic  Maps  of  Quality  Control  Sampling  Sites  of  the  Georgia 
Ecoregions  Proj  ect.  (cont. ) 


Ecoregion  75  -  Southern  Coastal  Plains 
Spatial  Quality  Control  Sampling  Sites 


Figure  45.  Spatial  Quality  Control  Sampling  Sites  for  Ecoregion  75. 
Ecoregion  75  -  Spatial  Quality  Control  Sampling  Site  List 


Ecoregion    Subecoregion        Site  ID 


Site  Name 


Impairment 
Status 


75e 


75 


75e  -  36 
75e  -  59 


Swain  Creek 
Ray  Branch 


Impaired 
Reference-GIS 


75f 


75f-44 


Springfield  Canal 


Impaired 


75j 


75J-41 


White  Branch 


Reference 
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Appendix  E.     Geographic  Maps  of  Quality  Control  Sampling  Sites  of  the  Georgia 
Ecoregions  Project,  (cont.) 

Ecoregion  75  -  Southern  Coastal  Plains 
Temporal  Quality  Control  Sampling  Sites 


Figure  46.  Temporal  Quality  Control  Sampling  Sites  for  Ecoregion  75. 


Ecoregion  75  - 

Temporal  Quality  Control  Sampling  Site  List 

Ecoregion    Subecoregior 

i      Site  ID 

Site  Name 

Impairment 
Status 

75h 

75h  -  70 

Pond  Fork 

Impaired 

75 

75J-2 
75j  -  24 

Trib.  to  Little  Ogeechee 

River 

Yellow  Bluff  Creek 

Impaired 
Impaired 
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